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Methods and Apparatus to Position a Mobile Receiver using Downlink Signals 

Field of the Invention: This invention relates to location finding and tracking of a 
Mobile Receiver (MR) using a Reference Receiver (RR) of known location. 

5 

Background of the Invention; 

Whenever a number of signals are avmlable for reception to the MR as well as 
to the RR, it is possible to locate the MR. Examples of radio signals that exist with 
worldwide coverage are satellite signals such as Global Positioning System (GPS), 

10 Iridium, Globalstar, etc. Examples of radio signals that can penetrate structures such as 
buildings are land-based systems such as cellular, PCS, paging, and ESMR, among 
others. Both satellite-based signals and terrestrial-based signals can be used to locate 
the MR by relying on their reception at both the MR and the RR. In this patent, we 
refer to such a location system as a handset-based Wireless Location System (WLS) (or 

15 Forward/downlink-based WLS), as opposed to a network-beised WLS (or 
Reverse/uplink-based WLS), which locates a Mobile Transmitter (MT) using a 
plurality of receivers that receive the signal transmitted from the MT. An advantage of 
a network-based WLS is that it can be passive (in that it can locate MR*s without 
modification or addition to wireless transmissions) and that no modifications to the 

20 hardware or operation of the MT are required for the WLS to locate the MT. On the 
other hand, a handset-based WLS has the advantage of requiring no modification of 
existing network infrastructure. The location in the handset-based WLS can be 
performed either at the MR, or at the RR, or at a Third Receiver (TR). When the 
location is performed at the MR, the RR is required to transmit the Reference 

25 Information to the MR. When the location is performed at the RR, the MR is required 
to transmit the Location Information to the RR. When the location is performed at the 
TR, both the MR and the RR are required to transmit their corresponding information 
to the TR. 

Radio signals can interfere with one another unless multiple access techniques 
30 are employed. Existing radio signals use any one of a number of multiple access 
techniques such as: Code Division Multiple Access (CDMA), Time Division Multiple 
Access (TDMA) and Frequency Division Multiple Access (FDMA). CDMA is used in 
IS-95, the North American standard for CDMA CTs, as well as in GPS. FDMA is used 
in AMPs, the North American standard for analogue CTs. TDMA is used in IS- 136, the 



North American standard for TDMA CTs, and in GSM, the European standard for 
TDMA CTs. Without loss of generality, we will use in this patent GPS signals as an 
example of a satellite-based radio signal and IS-9S signals as an example of a land- 
based radio signal. This choice does not preclude any other standard or any other 
5 multiple access technique. In this patent, we refer to a transmitter as a Base Station 
Transmitter (BST). It is important to note that a KR is not required in this patent if the 
BSTs are referenced to the same reference signal. 

The need for wireless location finding and tracking of MRs is ever increasing. 
Some of the possible services for locating MRs are: 

10 

1. Enhanced Emergency Services: In recent years pressure has been mounting for 
the development of technology to position Cellular Telephones (CTs). The 
primary driving force has been enhanced 911 (E911) services for wireless 
telephone subscribers. E911 services provide the 911 operator with information 

15 such as calling number, street address, and the primary subscriber's name. A 

ruling concerning E91 1 emergency calling systems was released by the FCC in 
October, 1996 (CC Docket No. 94-102) and clarified in December 1997. In the 
document, the FCC requires that by October 2001, a wireless 91 1 caller be located 
with a horizontal accuracy of 125 meters RMS (Root Mean Square). 

20 2. Tracking of Fraudulent Calls. 

3. Tracking of Stolen Vehiciesi Conventionally, to track stolen vehicles requires 
installing an RF tag such as in Bird, U.S. Patent No. 5,418,537 issued May, 23, 
1995, leaving it permanently on in each vehicle to be tracked and a new 
infrastructure for the BSs throughout the desired service area. An alternative 

25 solution is to either use existing CTs to take advantage of the existing cellular 

coverage and the broad availability of inexpensive CTs, or to use GPS receivers to 
take advantage of the availability and accuracy of satellite signals. 

4. Fleet Management for Courier and Transportation Businesses: Once again, to 
manage a fleet requires installing an RF tag such as in Song, U.S. Patent No. 

30 5,208,756 issued May, 4, 1993, and Sheffer et al., U.S. Patent No. 5,218,367 

issued June 8, 1993. Instead, existing CTs can be used to take advantage of the 
existing cellular coverage and the broad availability of inexpensive CTs. 
Alternatively, GPS receivers can be used to take advantage of the availability and 
accuracy of satellite signals. 



5. Location-sensitive billing. 

Without loss of generality, we will assume throughout this document that the 
intention of the handset-based WLS is to estimate the horizontal position of the MR as 
well as its horizontal velocity. In case both the vertical and the horizontal position of a 
5 MR are to be estimated, an extra independent measurement is required in addition to 
the minimum number required for horizontal positioning. 

Objects of the Invention : 

According to an aspect of the invention, there is provided a handset-based 
10 Wireless Location System (WLS) whereby an MR can be located using signals 
received by both the MR and an RR. More specifically, it is an object of the invention 
to estimate the static and kinematic positional information of a MR which receives 
signals transmitted from a number of BSTs. 

According to a further aspect of the invention, there is provided a system that 
15 corrects for the sources of error that affect the different location technologies for a 
handset-based WLS in a unique and novel way. 

According to a further aspect of the invention, there is provided a system that 
solves for the position of the MR using Least Squares where misclosures cind 
standardized residuals are used to flag the observations that might have a blunder. 
20 According to a further aspect of the invention, there is provided a system that 

uses Chaffee's method or Location On the Conic Axis (LOCA) to detect solution 
bifurcation. 

According to a further aspect of the invention, there is provided a system that 
uses LOCA and/or Plane Intersection to provide an initial position for Least Squares. 
25 According to a further aspect of the invention, there is provided a system that 

uses individually, or in combination, TDOA observations, hybrid TDOA observations 
and AOA observations to locate the MR. 

According to a further aspect of the invention, there is provided a system that 
resolves the AOA ambiguity at a cellular BS by using the TDOA information from two 
30 distant BSs or the AOA from several sectors at the same cellular BS. 

According to a further aspect of the invention, there is provided a system that 
uses the range information from two BSTs or the AOA from several sectors at the same 
cellular base station to resolve the AOA ambiguity at such a BST. 



Objects of the Invention : 

According to an aspect of the invention, there is provided a handset-based 
Wireless Location System (WLS) whereby an MR can be located using signals 
received by both the MR and an RR. More specifically, it is an object of the invention 
5 to estimate the static and kinematic positional information of a MR which receives 
signals transmitted from a number of BSTs. 

It is also an aspect of the invention to provide a system which corrects for the 
sources of error that affect the different location technologies for a handset-based WLS 
in a unique and novel way. More specifically, it is proposed to correct for frequency 
10 errors, clock offsets, overall group delays, multipath and interference. 

Yet another aspect of the invention is to estimate the Doppler shift (6fi) in each 
received signal at the MR in order to estimate the speed and Direction of Travel (DOT) 
of the MR using Frequency Difference Of Arrival (FDOA). This applies to 
modifications: I and IL 

15 Yet another aspect of the invention is to minimize the use of RRs in order to 

reduce the cost and complexity of a handset-based WLS. It is an object of the invention 
to produce a handset-based Wireless Location System (WLS) where an MR can be 
located without the need for an RR. This applies to all three modifications: 0, 1 and II. 

20 Objects of the Invention : 

According to an aspect of the invention, there is provided a handset-based 
Wireless Location System (WLS) whereby an MR can be located using signals 
received by both the MR and an RR. More specifically, it is an object of the invention 
to estimate the static and kinematic positional information of a MR which receives 

25 signals transmitted from a number of BSTs. 

Another contribution of this invention is the way a CT is requested to act as an 
MR in a handset-based WLS where a transmitter in a cellular BS is considered to be a 
BST. One should keep in mind that for a CT to act as an MR, a certain level of CT 
modification might be required. Since there is a large number of existing CTs 

30 worldwide, it is important to keep in mind that some levels of modification to the CT 
might preclude existing CTs (which is not a desirable outcome). We attempt to reduce 
the amount of CT modification for a desired WLS performance level. Three levels of 
CT modification have been identified in this patent: 
modification 0 which requires no modification to the CT; 
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modification I which requires only software modification to the existing CT; and 
modification II which requires both software and hardware modifications to the 
existing CT. 

It is an aspect of the invention to improve modification 0, 1 and II over existing 
5 methods. 

It is an aspect of the invention to provide a system that optimizes the threshold 
used in selecting the correlation peak, and to optimize the list of signals considered for 
correlation, for all modifications: 0, I and II. It is a ftirther aspect of the invention to 
optimize both the threshold as well as the list of signals considered for correlation 

10 based on minimizing HDOP. 

Yet another aspect of the invention is the use of the SR and Inverse SR 
algorithms over the correlation function (that is obtained by correlating the received 
signals at both the MR and RR with a clean version of themselves) in order to: 
resolve the correlation peaks beyond the Rayleigh resolution and to 

15 avoid using a threshold in selecting the correlation peak. 

According to another aspect of the invention, there is provided a system that 
optimizes the list of signals considered for correlation based once again on minimizing 
HDOP. This applies to modifications: I and II. 

Yet another aspect of the invention is a system for locating an MR by estimating 

20 the Phase Of Arrival (POA) of existing or generated tones contained in the received 
signals. The generated tones can be obtained by performing a nonlinear operation over 
the received signals at the MR. The phases of the existing or generated tones can be 
extracted using SR algorithms in order to reduce the effect of multipath by selecting 
the tone corresponding to the first arrival in time. This applies to modifications: I and 

25 II. 

Yet another aspect of the invention is use of the estimated POA to estimate 
either the TO A of the received signals or their AOA using Phase Difference Of Arrival 
(PDOA). This applies to modifications: I and 11. 

Yet another aspect of the invention is provision of a system which estimates the 
30 Doppler shift (5fi) in each received signal at the MR in order to estimate the speed and 
Direction of Travel (DOT) of the MR using Frequency Difference Of Arrival (FDOA). 
This applies to modifications: I and II. 

Yet another aspect of the invention is a system that minimizes the use of RRs in 
order to reduce the cost and complexity of a handset-based WLS. According a further 



aspect of the invention, there is provided a handset-based Wireless Location System 
(WLS) where an MR can be located without the need for an RR. This applies to all 
three modifications: 0, 1 and 11. 

According to a further aspect of the invention, there is provided a system that 
5 resolves the AOA ambiguity at a cellular BS by using the TDOA information from two 
distant BSs or the AOA from several sectors at the same cellular BS, 

According to a further aspect of the invention there is provided a system which 
uses the range information from two BSTs or the AOA from several sectors at the same 
cellular base station to resolve the AOA ambiguity at such a BST. 

10 

Objects of the Invention ; 

According to an aspect of the invention, there is provided a handset-based 
Wireless Location System (WLS) whereby an MR can be located using signals 
received by both the MR and an RR. More specifically, it is an object of the invention 
15 to estimate the static and kinematic positional information of a MR which receives 
signals transmitted from a number of BSTs. 

According to a further aspect of the invention, there is provided a system that 
corrects for the sources of error that affect the different location technologies for a 
handset-based WLS in a unique and novel way. 
20 According to a further aspect of the invention, there is provided a system that 

solves for the position of the MR using Least Squares where misclosures and 
standardized residuals are used to flag the observations that might have a blunder. 

According to a further aspect of the invention, there is provided a system that 
uses Chaffee's method or Location On the Conic Axis (LOCA) to detect solution 
25 bifurcation. 

According to a further aspect of the invention, there is provided a system that 
uses LOCA and/or Plane Intersection to provide an initial position for Least Squares. 

According to a further aspect of the invention, there is provided a system that 
uses individually, or in combination, TDOA observations, hybrid TDOA observations 
30 and AOA observations to locate the MR. 

According to a further aspect of the invention, there is provided a system that 
resolves the AOA ambiguity at a cellular BS by using the TDOA information from two 
distant BSs or the AOA from several sectors at the same cellular BS. 



wo 01/65271 PCT/CAOO/00224 



According to a further aspect of the invention, there is provided a system that 
uses the range information from two BSTs or the AOA from several sectors at the same 
cellular base station to resolve the AOA ambiguity at such a BST. 

5 Definitions: 

The following definitions take precedence over definitions for the same terms that can 
be found in the open literature. 

Mobile Receiver (MR); is a device, which can be portable or fixed, that can consist of 
a receiver alone or both a receiver and a transmitter. It can be a regular Cellular 

10 Telephone (CT), a PCS (Personal Communication Systems) telephone, a cordless 
telephone, a Personal Digital Assistant (PDA), a GPS receiver, or a combination 
thereof. It can be a radio tag or a wireless telephone that does not contain the audio 
portion of the telephone. It can also be a transmitter that transmits periodically over a 
given chaimel, or a receiver that receives Radio signals, or both. 

J 5 Reference Receiver CRR): is a device of known location with respect to the location of 
BSTs that are transmitting radio signals intended for reception by both the MR and the 
RR. An RR can be portable or fixed. It can consist of a receiver alone or both a receiver 
and a transmitter. The RR estimates the TOA, POA, SOA and/or FOA of signals 
transmitted by the BSTs. It either transmits the Reference Information to the MR or to 

20 the TR, or it receives the Location Information from the MR. In some unique situations 
when the BSTs are referenced to a common signal, the RR is not required. 
Reference Information; is the estimated TOA, POA, SOA and/or FOA of the signals 
transmitted by a number of BSTs at the RR as well as the location of such BSTs. 
Location Information; is the estimated TOA, POA, SOA and/or FOA of the signals 

25 transmitted by a number of BSTs at the MR. 

Third Receiver (TR>; is a device that receives Location Information from the MR as 
well as Reference Information from the RR in order to estimate the location of the MR. 
It can be a receiver only, or a transmitter and a receiver. The Mobile Switching Center 
(MSC) can play the role of the TR. 

30 Base Station Transmitter (BST); is a device of known location with respect to other 
BSTs at a given time. A BST can be portable or fixed. It can consist of a transmitter 
alone or both a receiver and a transmitter. It can be a regular cellular Base Station (BS), 
a regular satellite transceiver, a GPS satellite, a Loran-C transmitter, a dedicated 
transmitter, a transponder, or any other type of transmitter/transceiver combination. 
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8 

Host: is a processor that processes the Location Information together with the 
Reference Information to estimate the location of the MR. The function of the Host can 
be played by the MR, the RR or a TR. The host can also be responsible for location 
services such as fleet management, location sensitive billing, etc. 
5 Network-based WLS rFonvard/Downlink-based WLS): locates a Mobile Receiver 
(MR) using a plurality of BSTs that transmit signals received by both the MR and the 
RR. 

Handset-based WLS (Reverse/uplink-based WLS): locates a Mobile Transmitter 
(MT) using a plurality of receivers that receive the signal transmitted from the MT. 

10 Cell: is a geographical area serviced by a cellular Base Station (BS). 

Sectorized cell: is a cell that is made of several non-overlapping sectors. Each sector 
can be considered as an independent cell to be serviced by an independent cellular BS. 
However, sectors in a cell are usually serviced by the same cellular BS in order to 
minimize cost and complexity. We refer to such a BS as a sectorized BS. 

15 Pilot Channel: is the channel to be used by each sector to broadcast a (locally unique) 
pilot signal. 

Pilot Signal: is defined as a sequence known to the MR. 

As defined in the IS-95 standard, this sequence is a code repeated exactly 75 times 
every 2 seconds, and is known as the "pilot FN sequence'' or the "short code". 

20 Although all sectors broadcast the same pilot signal, the transmissions from each base 
station have different synchronization or phases, which distinguish the pilot signals. 
Phase increments are specified as multiples of 64 chips, where 1 chip is 1/1.2288 
microseconds. The offset in units of 64 chips, is referred to as the pilot PN offset. 
Pilot increments in a network are typically 3 or 4 PN offsets. A PN offset of 4 

25 corresponds to 256 chips, or 208.3 microseconds. 

As defined in the AMPs standard or in the IS-136/GSM standard, this sequence is a 
Digital Colour Code (DCC) that is repeated at specific intervals of time. 
Super-Resolution (SR) Algorithm; is an operation that transforms a time domain 
signal, s(t), to a frequency domain signal, S(f), in such a way that the frequency domain 

30 signal, S(f), has a better resolution than the resolution offered by the Fourier transform, 
i.e. S(f) has a better resolution than the Fourier transform F{s(t)}, of s(t). Vice- versa, a 
SR algorithm transforms a frequency domain signal, S(f), to a time domain signal, s(t), 
with a better resolution than the resolution offered by the inverse Fourier transform, i.e. 
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s(t) has a better resolution than the inverse Fourier transform F-^{S(f)}, of S(f)). 
Examples of SR algorithms are well known in the literature and include: 
MUSIC/Root MUSIC, ESPRIT, Auto Regressive Moving Average, Minimum 
Variance, MUSIC using Higher Order Statistics, ESPRIT using Higher Order 
5 Statistics, Auto Regressive Moving Average using Higher Order Statistics, or 
Minimum Variance using Higher Order Statistics. 

Frequency-domain Ravleigh Resolution; is the frequency-domain resolution offered 
by the Fourier transform. 

Time-domain Ravleigh Resolution: is the time-domain resolution offered by the 

10 inverse Fourier transform. 

Inverse SR algorithm; is an algorithm which processes a time domain signal in order 
to improve its time-domain resolution over the conventional time-domain Rayleigh 
resolution. A preferred embodiment of the inverse SR algorithm comprises 
a conventional time domain correlator, a time domain window, a Fourier Transform, 

15 a frequency domain window, a frequency domain equalizer, or a processor which 
performs a SR algorithm in order to resolve the TOAs in the received signal at a given 
MR. 

Altematively, an inverse SR algorithm can process a frequency domain signal in 
order to improve its frequency-domain resolution over the conventional frequency- 

20 domain Rayleigh resolution. In this case, a preferred embodiment of the inverse SR 
algorithm comprises a frequency domain correlator, a frequency domain window, an 
inverse Fourier Transform, a time domain window, a time domain equalizer, or 
a processor which performs a SR algorithm in order to resolve the FOAs in the 
received signal at a given MR. 

25 Sector Antennas; are cellular antennas that exist in some BS that service more than 
one sector. Each sector has a dedicated antenna (or set of antennas if diversity is 
applied). 

Fourier Transform-Based Filters; are filters that: Fourier Transform the time domain 
signal, then window the transformed signal over a given band, and 
30 Inverse Fourier Transform the windowed signal. 

Description of an IS-95 CT receiving Pilot signals ; 

The CT measures the quality of received pilot signals by computing the ratio between 
the received energy of a pilot (Ec) to the total received power by the mobile (Iq). This 



10 



measurement is referred to as the EJlo of the pilot. This measurement is closely related 
to the signal strength of the pilot. An IS-95 CT uses the Ec/Io measurement of a pilot to 
estimate the suitability of the communications link between itself (the CT) and that 
particular sector in a BS. For example, if a CT is in a mode in which it may only 

5 commimicate with one sector at a time, it may choose to communicate with the sector 
which has the highest Ec/Io measured by the CT. 

An IS-95 CT may include a searcher element and multiple demodulator fingers 
as shown in Figure 12. A searcher is used to search for pilot signals by measuring the 
Ec/Io of a pilot signal at different delays, A demodulator finger tracks a pilot signal, 

10 and may demodulate an information signal associated with said pilot signal. 

In addition to transmitting a pilot signal, a base station may also transmit one or 
more information signals associated with the pilot. These information-bearing signals 
are transmitted at the same radio frequency as their associated pilot signal and are 
synchronized with the pilot signal. Thus, an information signal will undergo 

15 propagation effects (attenuation, phase shift, time delays, multipath, frequency shift) 
nearly identical to those of its associated pilot signal. A mobile then tracks the 
propagation effects of a pilot to aid in demodulating information from an associated 
signal, or code channel. As well, the mobile can estimate the reliability of information 
demodulated from a code channel from the Ec/I© quality measurement of the associated 

20 pilot signal. 

In general, a CT searches for pilot signals using one or more searchers, and then 
based on the searcher results, assigns demodulator fingers to demodulate information 
from code channels associated with suitable pilots. This process is ongoing so that the 
finger assignment adapts to the changing conditions typical of a CT cellular 
25 environment. 

During normal operation, the searchers and demodulator fingers in an IS-95 CT 
are usually assigned according to algorithms which optimize the communications 
performance of the phone, and the network as a whole. These assignment algorithms 
will limit the performance of estimation of the CT location, as they will often ignore 
30 pilot signals usefiil for the location estimation. We prefer to force the CT to enter into a 
state optimized for location. 

Summary of the Invention 
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11 

In order to accomplish the aforementioned objects and intentions, there are 
provided inventive contributions which are claimed in this patent document. 

Therefore, According to an aspect of the invention, there is provided a method 
5 of estimating the location of a Mobile Receiver (MR) comprising the steps of: 

receiving a plurality of signals from a plurality of Base Station Transmitters 
(BSTs) by the MR; and 

performing a nonlinear operation on the plurality of received signals at the MR 
in order to generate at least one tone per received signal. 
10 Additional refinements of this aspect of the invention are set forth in claims 2-5 

attached hereto, which are hereby incorporated in this summary by reference. 

According to an aspect of the invention, there is provided a method of 
estimating the location of a Mobile Receiver (MR) comprising the steps of: 

receiving a plurality of signals from a plurality of Base Station Transmitters 
1 5 (BSTs) by the MR; and 

estimating the location of the MR by estimating the phases of existing tones in 
each received signal. 

Additional refinements of these aspects of the invention are set forth in claims 7- 
17 attached hereto, which are hereby incorporated in this summary by reference. 
20 According to an aspect of the invention, there is provided a method of estimating 

the location of a MR from a plurality of signals transmitted from a plurality of BSTs, 
where one or more antennas is used at the MR and in which location of the MR is 
achieved using a technique selected from the group consisting of:: 
AOA positioning (multi-angulation), 
25 TOA positioning (circular multilateration). 

Range positioning (circular multilateration), 
TDOA positioning (hyperbolic multilateration), 
TDOA positioning (circular multilateration with TDOAs), 
AOA/TDOA positioning (muhi-angulation/hyperbolic 

30 multilateration). 
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AOA/TDOA positioning (multi-angulation/circular multiiateration 
with TDOAs), 

AOA/Range positioning (multi-angulation/circular multiiateration), 
AOA/Range positioning (multi-angulation/circular multiiateration), 
and 

combining the signals from all antennas for location purposes. 
Additional refinements of this aspect of the invention are set forth in 



claims 1 9, 20 and 22 attached hereto, which are hereby incorporated in this summary by 
reference 

According to an aspect of the invention, there is provided a method of 
estimating the location of a MR from a plurality of signals transmitted from a plurality 
of BSTs, using equation (29) and equation (30) (defined below). 

According to an aspect of the invention, there is provided a method of 
estimating the location of a MR wherein location of the MR is achieved based on a 
solution selected from the group consisting of: 

solution I, 

solution II, 

solution III, and 

solution IV (as described below). 

According to an aspect of the invention, there is provided a system to estimate 
the location of a Mobile Receiver (MR) comprising of: 

a receiver at the MR for receiving a plurality of signals from a plurality of Base 
Station Transmitters (BSTs); and 

a computing means coimected to the receiver by a communication linkfor 
performing a nonlinear operation on the plurality of received signals in order to 
generate at least one tone per received signal. 

Additional refinements of this aspect of the invention are set forth in claims 26 
and 27 attached hereto, which are hereby incorporated in this summary by reference 
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According to an aspect of the invention, there is provided a system to estimate 
the location of a Mobile Receiver (MR) comprising of: 

a receiver at the MR for receiving a plurality of signals from a plurality of Base 
Station Transmitters (BSTs); and 

a computing means connected by a communication link to the receiver for 
estimating the location the MR by estimating the phases of existing tones in each 
received signal. 

Additional refinements of this aspect of the invention are set forth in claims 2-5 
attached hereto, which are hereby incorporated in this siunmary by reference 

According to an aspect of the invention, there is provided a system to estimate 
the location of a MR from a plurality of signals transmitted from a plurality of BSTs, 
where one or more antennas are used at the MR and wherein the system comprises 
means for implementing a technique selected from the group consisting of: 

AOA positioning (multi-angulation), 

TOA positioning (circular multilateration), 

Range positioning (circular multilateration), 

TDOA positioning (hyperbolic multilateration), 

TDOA positioning (circular multilateration with TDOAs), 

TOA/TDOA positioning (multi-angulation/hyperbolic 

multilateration), 

AO A/TDO A positioning (multi-angulation/circular multi lateration 
with TDOAs), 

AOA/Range positioning (multi-angulation/circular multilateration), 
AOA/Range positioning (multi-angulation/circular multilateration), 
and 

combining the signals from all antermas for location purposes. 
Additional refinement of this aspect of the invention is set forth in claim 29 and 
31 attached hereto, which are hereby incorporated in this sununary by reference. 
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According to an aspect of the invention, there is provided a system to estimate 
the location of a MR from a plurality of signals transmitted from a plurality of BSTs, 
using equation (29) and equation (30). 

According to an aspect of the invention, there is provided a system to estimate 
5 the location of a MR from a plurality of signals transmitted from a plurality of BSTs, 
based on any one of the following solutions: 

solution I, 

solution II, 

solution III, or 
10 solution IV. 

According to an aspect of the invention, there is provided a system to estimate 
the location of a Mobile Receiver (MR), the system comprising a Receiver that 
receives a plurality of signals transmitted from a plurality of BSTs; 



Additional refinements of this aspect of the invention are set forth in claims 34- 
44 and 49-54 attached hereto, which are hereby incorporated in this summary by 
reference. 

According to an aspect of the invention, there is provided a system to estimate 
25 the location of a MR using a Reference Receiver havmg a LO, the system comprising: 

an estimator for estimating the carrier offset (due to Doppler, 8fi, and due to LO 
offset, Afi) of a plurality of signals transmitted by a plurality of BSTs and received by 
the RR. 



20 



15 



wherein the Receiver includes an IF-sampling Receiver comprising: 
a number of IF stages to convert a received RF signal ri(t) to an analog IF 
signal of IF frequency fi; 
a pre-aliasing filter; 

an Analog-to-Digital (A/D) converter that converts the IF analog signal to an IF 
digital signal; and 

a processor for processing the IF digital signal. 
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Additional refinements of this aspect of the invention are set forth in claims 45- 
54 and 60 attached hereto, which are hereby incorporated in this summary by reference. 

According to an aspect of the invention, there is provided a system to estimate 
the location of a MR having a LO, the system comprising: 
5 an estimator for estimating the carrier offset (due to Doppler, 5fi, and due to LO 

offset, Afi) of a plurality of signals transmitted by a plurality of BSTs and received by 
. the MR. 

Additional refinements of this aspect of the invention are set forth in 
claims 56-59 and 61-66 attached hereto, which are hereby incorporated in this summary 
10 by reference 

According to an aspect of the invention, there is provided a method of 
estimating the location of a MR, the method comprising the steps of: 

receiving a plurality of signals transmitted from a plurality of BSTs at a 
Receiver; 

15 converting the received RF signal ri(t) to an IF signal of IF frequency fi in an 

IF-sampling Receiver; 

filtering the IF signal with a pre-aliasing filter; 

converting the filtered IF analog signal to an IF digital signal in an Analog-to- 
Digital converter; and 
20 processing the IF digital signal in a signal processor. 

Additional refinement of this aspect of the invention is set forth in claim 68 
attached hereto, which are hereby incorporated in this summary by reference. 

According to an aspect of the invention, there is provided a method of 
estimating the location of a MR, the method comprising the steps of: 
25 estimating the carrier offset (due to Doppler, Dfj, and due to LO offset, Dfi) of a 

plurality of signals transmitted by a plurality of BSTs and received by a RR. 

According to an aspect of the invention, there is provided a method of 
estimating the location of a MR, the method comprising the steps of: 
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estimating the carrier offset (due to Doppler, 5fi, and due to LO offset, Afi) of a 
plurality of signals transmitted by a plurality of BSTs and received by the MR. 

According to an aspect of the invention, there is provided a method of 
estimating the location of a MR from a plurality of signals transmitted from a plurality 
5 of BSTs, the method comprising the steps of: 

characterizing the relative group delay, Atj - Atj, at a RR between a signal 
received from the i^** EST and a signal received from the j'*^ EST prior to operation; 

removing the characterized relative group delay, Ati - Atj, from the estimated 
Time Difference Of Arrival, x\ - xj, at the MR between a signal received from the i^^ 
10 EST and a signal received from the EST to yield a corrected Time Difference of 
Arrival; and 

finding the location of the MR using the corrected Time Difference of Arrival. 

Additional refinements of this aspect of the invention are set forth in claims 72- 
74 attached hereto, which are hereby incorporated in this summary by reference. 
15 According to an aspect of the invention, there is provided a method of 

estimating the location of a MR in location state 0 or location state 1 from a plurality of 
signals transmitted from a plurality of BSTs, wherein the transfer of the positional 
information from the Host to a user is achieved over the internet. 

Additional refinements of this aspect of the invention are set forth in claims 76- 
20 78 attached hereto, which are hereby incorporated in this summary by reference. 

According to an aspect of the invention, there is provided a method of 
estimating the location of a MR from a plurality of signals transmitted from a plurality 
of BSTs, comprising: 

using Least Squares to solve for the positional information of the MR and 
25 using TDOA observations from distant BSTs to resolve for ambiguities in an AOA 
solution; 

using AOA observations from different sectors at the same cellular BS to 
resolve for ambiguities in an AOA solution; and 
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using range observations from different BSTs to resolve for ambiguities in an 
AOA solution. 

Additional refinement of this aspect of the invention is set forth in claim 80 
attached hereto, which are hereby incorporated in this summary by reference 

According to an aspect of the invention, there is provided a system for 
estimating the location of a MR from a plurality of signals transmitted from a plurality 
of BSTs, the system including a host, and the system comprising; 

the MR being adapted to receive some or all the following information from a 
remote source and to apply said information to the execution of a location state: 

identification of base station pilot signals for the mobile to search for pilot 
detection threshold; 

duration of time for searcher to accumulate pilot energy at each delay offset; 

search window time length or duration; 

spacing of successive sampling points within the search window; 
frequency at which to conduct the search for pilot signals; 
cellular system which the mobile should acquire; 

the number of times the mobile should repeat all or part of the location 
procedure; 

correlation interval; 
locations of base stations; and 

the MR being adapted to transmit one or more of the following information to 
the Host: 

times of arrival of pilot signals; 

identification of the pilot signals reported by said device; 
received signal qualities of pilot signals; 
pilot search results; and 

carrier frequencies of reported pilot signals. 



18 



According to an aspect of the invention, there is provided a system for 
estimating the location of a MR from a plurality of signals transmitted from a plurality 
of BSTs, the system comprising: 

first computing means for characterizing the relative group delay, Atj - Atj, at a 
5 RR between a signal received from the i^ BST and a signal received from the j^*" BST 
prior to operation; and 

second computing means for removing the characterized relative group delay, 
At| - Atj, from the the estimated Time Difference Of Arrival, x\ - xj, at the MR between a 
signal received from the i^*^ BST and a signal received from the BST. 
10 According to an aspect of the invention, there is provided a system 

for estimating the location of a MR in location state 0 or in location state 1 
from a plurality of signals transmitted from a plurality of BSTs, wherein the 
transfer of the positional information from the Host to a user is achieved over 
the internet. 

15 Additional refinement of this aspect of the invention is set forth in claim 84 

attached hereto, which are hereby incorporated in this summary by reference. 

According to an aspect of the invention, there is provided a system for 
estimating the location of a MR from a plurality of signals transmitted from a plurality 
of BSTs, where: 

20 Least Squares is used to solve for the positional information of the MR; 

TDOA observations from distant BSTs are used to resolve for ambiguities in an 
AOA solution; 

AOA observations from different sectors at the same cellular BS are used to 
resolve for ambiguities in an AOA solution; and 
25 AOA observations from different sectors at the same cellular BS are used to 

resolve for ambiguities in an AOA solution. 

Additional refinement of this aspect of the invention is set forth in claim 86 
attached hereto, which are hereby incorporated in this summary by reference. 
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According to an aspect of the invention, there is provided a method to estimate 
the location of a MR from a plurality of signals transmitted from a plurality of BSTs; 

the MR receiving some or all the following information from a remote source 
and applies said information to the execution of a location state: 

identification of base station pilot signals for the mobile to search for pilot 
detection threshold; 

duration of time for searcher to accumulate pilot energy at each delay offset; 

search window time length or duration; 

spacing of successive sampling points within the search window; 
frequency at which to conduct the search for pilot signals; 
cellular system which the mobile should acquire; 

the number of times the mobile should repeat all or part of the location 
procedure; 

correlation interval; 
locations of base stations; and 

wherein some or all of the following information is transmitted from the MR to 
the Host: 

times of arrival of pilot signals; 

identification of the pilot signals reported by said device; 
received signal qualities of pilot signals; 
pilot search results; 

carrier frequencies of reported pilot signals. 

According to an aspect of the invention, there is provided a method of 
estimating the location of a MR from a plurality of signals transmitted from a plurality 
of BSTs, where misclosures and/or standardized residuals are used to flag observations 
that might contain a blimder. 

Additional refinements of this aspect of the invention are set forth in 
claims 89-97, 106-112 and 114 attached hereto, which are hereby incorporated in this 
sununary by reference. 
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According to an aspect of the invention, there is provided a method of 
estimating the location of a MR from a plurality of signals transmitted from a plurality 
of BSTs, where a combination of AOA and TDOA observations are used such that the 
MR can be located with as few as two AOA independent observations and one TDOA 
5 observation. 

Additional refinements of this aspect of the invention are set forth in claims 99 
attached hereto, which are hereby incorporated in this summary by reference. 

According to an aspect of the invention, there is provided a method of 
estimating the location of a MR from a plurality of signals transmitted from a plurality 
10 of BSTs, wherein a combination of AOA and TOA observations are used such that the 
MR can be located with as little as one AOA observation and one TOA observation. 

Additional refinements of this aspect of the invention are set forth in claims 101- 
102 attached hereto, which are hereby incorporated in this summary by reference. 

According to an aspect of the invention, there is provided a method of 
15 estimating the location of a MR from a plurality of signals transmitted from a plurality 
of BSTs, where a combination of TOA observations are used such that the MR can be 
located with as little as two independent TOA observations. 

Additional refinements of this aspect of the invention are set forth in claims 104- 
1 05 attached hereto, which are hereby incorporated in this summary by reference 
20 According to an aspect of the invention, there is provided a method of 

estimating the location of a MR from a plurality of signals transmitted from a plurality 
of BSTs, which quantifies the effect of geometry on speed and direction of travel 
estimation using the design matrix in equation (51). 

According to an aspect of the invention, there is provided a system to estimate 
25 the location of a MR from a plurality of signals transmitted from a plurality of BSTs, 
where misclosures and/or standardized residuals are used to flag observations that 
might contain a blunder. 
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Additional refinements of this aspect of the invention are set forth in 
claims 116-118 and 122-125 attached hereto, which are hereby incorporated in 
this summary by reference. 

According to an aspect of the invention, there is provided a system to estimate 
5 the location of a MR from a plurality of signals transmitted from a plurality of BSTs, 
where a combination of AOA and TDOA observations are used such that the MR can 
be located with as few as two AOA independent observations and one TDOA 
observation. 

According to an aspect of the invention, there is provided a system to estimate * 
10 the location of a MR from a plurality of signals transmitted from a plurality of BSTs, 
where a combination of AOA and TOA observations are used such that the MR can be 
located with as little as one AOA observation and one TOA observation. 

According to an aspect of the invention, there is provided a system to estimate 
the location of a MR from a plurality of signals transmitted from a plurality of BSTs, 
15 where a combination of TOA observations are used such that the MR can be located 
with as little as two independent TOA observations. 

Additional refinements of this aspect of the invention are set forth in claims 122- 
123 attached hereto, which are hereby incorporated in this summary by reference. 

According to an aspect of the invention, there is provided a system to estimate 
20 the location of a MR from a plurality of signals transmitted from a plurality of BSTs, 
which quantifies the effect of geometry on position estimation, with both AOA and 
TDOA observations, using the design matrix in equation (46). 

According to an aspect of the invention, there is provided a system to estimate 
the location of a MR fi*om a plurality of signals transmitted fi-om a plurality of BSTs, 
25 which quantifies the effect of geometry on speed and direction of travel estimation 
using the design matrix in equation (51). 

Description of the Figures ; 

Preferred embodiments of the invention will be described with reference to the 
30 Figures, by way of example, and without intending to limit the invention to the 
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specific embodiments disclosed, in which Figures like reference characters denote 
like elements, and in which: 

Figure 1 illustrates the broadcast of a signal, Si(t) (103), by the i^ BST (102) 
where To,i is the Time of Transmission (TOT) of the signal. A MR (101) receives the 
5 signal, together with other signals transmitted from a number of BSTs. 

When the BST (102) corresponds to an IS-95-based sector in a BS, the signal, 
Si(t) (103), can be a pilot signal, a sync signal, a paging signal or a Forward Traffic 
signal. When Sj(t) is a pilot signal, it consists of a single code, repeated exactly 75 times 
every 2 seconds, with a period of 32,768 chips. The pilot code is known as the "pilot 

10 PN sequence'' or the "short code". Although in IS-95 all sectors broadcast the same 
short code, the pilot transmitted by each sector is time-offset by an integral multiple of 
64 chips from any other pilot signal which could forseeably be simultaneously received 
by a MR. As the code has a period of 32,768 chips, and pilots are spaced by 64 chips, 
there are 512 unique pilot phases possible. 

15 When the BST (102) corresponds to a GPS-based satellite, the signal, Si(t) 

(103), can be either a Coarse Acquisition (CA) signal or a Precision signal. When Si(t) 
is a CA signal, it consists of a unique Gold code per satellite (i.e. up to 24 distinct 
codes), repeated exactly 1000 times every 1 second, with a period of 1023 chips, with 
BPSK overhead messages at 50 bits/second. The codes are known as the "CA-codes". 

20 When the BST (102) corresponds to an AMPs-based sector or an IS-136/GSM- 

based sector in a BS, the signal, Si(t) (103), can be a Forward Control Channel (FOCC) 
signal or a Forward Voice Channel (FVC) signal. When Si(t) is a FOCC signal, it 
consists of Digital Colour Codes (DCCs) with overhead messages regarding the 
channel numbers of the paging channels. Each DCC consists of known data, unique to 

25 its corresponding servicing sector, interspersed with overhead messages. 

Following Figure 1, Figure 2 illustrates the broadcast of the signal Si(t) by the i^^ 
BST (201) at Time 'xo,i' and its reception by the MR (202) at Time of Arrival: 'xi*, 

for i=l, 2, 3. 

When the BST is IS-95-based or GPS-based, the time To,i is the time of 
30 transmission of a given pilot PN sequence while xi is the time of reception of such a 
PN sequence. When the BST is AMPs-based, the time xo,i is the time of transmission 
of a given DCC sequence while xj is the time of reception of such a sequence. In order 
to solve for the horizontal coordinates (x,y) of the MR (202), the minimum number of 
BSTs required depends on the positioning technology used. 
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Figure 3 illustrates Angle of Arrival for several wavefronts from BSTs at a MR, 
and assumes that the MR is far from BSTj and BSTj with respect to the baseline 
between the two BSTs. Such an assumption implies that the transmitted wavefronts are 
planar, i.e. yi = Yj = yij, where: 
5 • Yi is the Angle Of Arrival (AO A), of signal Si(t) from the i**" BST to the MR, 

• Yj is the AOA of signal Sj(t) from the j^'' BST to the MR, and 

• YiJ is defined as the angle formed between: 

1 . the line joining the i* BST to the j*^ BST (which we refer to as Hney); 
and, 

10 2. the line joining the MR (301) and the middle of lineij between the two 

BSTs (which we refer to as linccij) 
in a clockwise manner from lineo^j to lineij. 

Figure 4 illustrates the Direction of Travel (DOT) *<p' (401) relative to Northing (in a 
15 clockwise mamier from Northing) and the speed v (402) of the MR (404) relative to 
the BSTs (403). '<p' and v represent the velocity v of the MR with respect to the BSTs 
(403). In a cellular system, the BSTs are stationary which is not the case in a satellite- 
based system. 

Without loss of generality. Figure 5 illustrates the two-dimensional (horizontal) 
20 Locus of Position (50 1 ) for TDOA2, 1 which is defined as 

TDOA2,, = (x2-x,2)-(T,-V,) 

= - V(X2 - X)' 4- (y^ - y)' - - V(xi - x)^ + (y . - y)^ 
c c 

(1) 

25 where c is the speed of propagation, (Xj , y| , Zj ) are the coordinates of the antenna of 
BSTi, (x2,y2jZ2) are the coordinates of the antenna at BST2, and (x,y) are the 
coordinates of the MR (502). This is achieved using TDOA-based Hyperbolic 
Multilateration (as shown in Turin, G. L. et al., "A Statistical Model of Urban 
Multipath Propagation," IEEE Transactions on Vehicular Technology, Vol. VT-21, 

30 No. 1, February 1972, and as shown in Smith, J.O. et al., "Closed-Form Least-Squares 
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Source Location Estimation from Range-Difference Measurements," IEEE 
Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP-35, No. 12, 
December, 1987). 

In the case where the three-dimensional coordinates of the MR are required, we 
5 have to solve for (x,y,2) using 



Without loss of generality. Figure 6 illustrates the two intersecting locii of Position for 
10 the two-dimensional coordinates (x,y) of the MR (603) based on TDOA2,l (601) and 
TDOA3J (602). In other words, it is possible to solve for (x,y) as the intersection 
between the two trajectories obtained using three Times of Arrival (TOAs) (after 
choosing the correct side of each hyperbola). In order to solve for (x,y,z) we require 
four TOAs. 

15 Figure 7a illustrates the Second Stage of the WLS where the signal, Si(t), is 

received by the MR (701). The MR (701) processes the received signal, ri(t), to 
estimate its own coordinates (x,y), speed v and DOT <p with the assistance of a 
Reference Receiver (RR) (702) of known location. The function of the RR is to 
estimate the Reference Information regarding the BSTs and to transmit it to the MR 

20 (701). We will refer to such a system as Option I. 

Figure 7b illustrates an alternative of the Second Stage of the Wireless location 
system where the signal, Si(t), is received by the MR (703). The MR processes the 
received signal, ri(t), and transmits the Location information to the RR (704). In turn, 
the RR estimates the coordinates (x,y) of the MR, its speed v and DOT <p, based on the 

25 Location Information transmitted by the MR (703) and based on the Reference 
Information estimated by the RR (704). We will refer to such a system as Option II. 

Figure 7c illustrates yet another alternative of the Second Stage of the Wireless 
location system where the signal, Si(t), is received by the MR (705). The MR processes 
the received signal, ri(t), and transmits the Location Informatioii to a Third Receiver 



TDOA2., = (X2 - 2 ) - (X| - x^^, ) 
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(TR) (706). In turn, the RR (707) estimates the Reference Information regarding the 
BSTs and transmits it to the same Third Receiver (TR) (706). The TR estimates the 
coordinates (x,y) of the MR, its speed v and DOT <p, based on the Location Information 
transmitted by the MR (705) and based on the Reference Information transmitted by 

5 the RR (707). We will refer to such a system as Option III. 

Figure 8 illustrates the description of Design I for the MR. The received RF 
signal ri(t) transmitted from the i^^ BST is initially filtered by an RF Band Pass Filter 
(BPF) (802), amplified by an RF amplifier (803), down-converted by one or more 
Intermediate Frequency (IF) stages (810). Each IF stage consists of a mixer (804) with 

10 an LO (808) generated using an RF synthesizer (809), an IF BPF (805) and an IF 
amplifier (806). After the IF stages, the resulting analog IF signal (807) of IF frequency 
fi is filtered by a pre-aliasing filter (811), sampled by an A/D (812) and filtered again 
using digital filters (813) in order to generate a digital IF signal (814) of center 
frequency f2, where f2 < fj. An optional Direct Digital Converter (DDC) (815) is used 

15 to down-convert the digital IF signal (814) to Baseband. The resulting signal is then 
processed by a Digital Signal Processor (DSP) (817) in order to estimate the Frequency 
Of Arrival (FOA) of the received signal, ri(t). A fiinction of the estimated FOA is fed- 
back to a Direct Digital Synthesizer (DDS) (819) which controls the RF synthesizer 
(809). The DSP (817) also provides a reference signal (820) to the DDS (819). A 

20 preferred source for the reference signal (820) is one that is common to all BSTs such 
as a GPS signal. 

Figures 9a, 9b, 9c and 9d illustrate a flow diagram for the positioning strategy. 
The positioning strategy consists of a number of algorithms (Least Squares, Location 
On the Conic Axis (LOCA), Plane Intersection) and numerous decisions. Given a set 
25 of observations, there are 4 possible outcomes: two LS position solutions (908), two 
LOCA position solutions (909), one LS position solution (923, 936), and no position 
solution (927). 

Figure 10 presents the preferred embodiment of an algorithm by which a MR 
can obtain TOA estimates of multiple IS-95 pilot signals. A search is made for each 
30 pilot over multiple time lags. If a pilot is detected, then it is fiirther observed to obtain 
a TOA estimate of refined accuracy. 

Figure 1 1 shows correlation values which might be obtained when searching for 
a particular pilot over multiple time lags, or delay offsets. Ec/Io is the ratio of average 
pilot chip energy to the total received power. A higher value of Ec/Io indicates greater 
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likelihood of arrival of the pilot signal at a particular delay offset. In this figure, the 
delay offsets are expressed in IS-95 chips (approximately 0.814 |j,s). In this figure, the 
horizontal dotted line indicates a threshold (-15 dB in this case) for detection of the 
pilot arrival. An Ec/Io value greater than this threshold is considered to indicate an 

S arrival of the pilot signal at the MR, at that particular time lag. 

Figure 12 indicates a possible architecture for the embodiment of the MR as a 
CDMA CT, which may be used to implement the algorithm of Figure 10. Under the 
direction of the Controller, the Searcher searches for pilot signals at different time lags 
of the received signal of the MR. The Demodulator Fingers are assigned to different 

10 pilot signals. When assigned to a pilot, a Demodulator Finger will slew its timing to 
track the pilot, and may demodulate data transmitted synchronously with that pilot. 
This architecture is presented in U.S. Patent number 5,764,687 Mobile Demodulator 
Architecture for a Spread Spectrum Multiple Access Communications System. 
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DETAILED DESCRIPTION OF PREFERJRED EMBODIMENT.^ OP 
THE INVENTION 

Theory : 

In a WLS, many factors affect the system perforaiance: 

1. RF shadowing and flat fading, 

2. jfrequency offsets (including LOs drift and Doppler Shifts), 

3. clock errors, 

4. time delays, 

5. noise» 

6^ multipath (selective feding), 

7. interference; 

8. geographical geometry of the BSTs relative to the intended MR. 

Each factor degrades the estimated location of the MR depending on the 
technolo^ employed for extracting the independent equations required for 
location. 

KE Transmission; More specifically, the Low Pass (LP) equivalent signal, 
^ (t) , traiisjxiitted from flie i.^ BST can be modeled as 

'^(t) = e"^'^<''*^>'-^ ei<^«Af;**«'> p,(t--T,^.) (3) 
and the.RF transmitted signal, Si(t), (see Figure 1) can be expressed as 

Si(t) = Re{e^^<^^(t)} (4) 

where 

• Re { . } denotes a real part operation; 

• ^ is the carrier frequency, 

• Afi is the frequency offset from f^ (usually unknown to the MR since it 
depends on the BST; it is to be estimated by the RR), 

• i = Time Of Transmission of signal p^t) (also unknown to the MR), 

• Q is the phase of the frequency f^ + Afj (also unknown to the MR), 

• j = V^,and 

• Pi(t) is a band-limited baseband signal transmitted by the i'** BST. 
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RF Reception: Then, the LP equivalent signal, ^(t). received by the MR and 
corresponding to ^(t) can be modeled as 

T(t) = G R:" e^^^'^*^^ ej'"^*''*'^'^' pjCt-tj + T„i - Ati) + Wi(t) + MPi(t)+ ^(t) 

(5) 

and the RF received signal, rjCt), (see Figure 2) can be expressed as 

ii(t) = Re{e^''^''^(t)i (6) 

where 

Re { . } denotes a real part operatioi^ 
i Gi represents the gain (real) due to the transmitting and receiving antennas: 
from the i* BST to the MR antenna (i.e. a fimction of the two antenna 
patterns), 

Rp represents the attenuation (real) due to the propagation channel, where 
n is a real number (usually between 2 and 4) that depends on the 
channel and 

, Rj is the. range bctwseeh the i*^ BSt and the MR, 



(7) 



is the phase of the received RF sigiial at the MR, 
where: 

Sf; is the frequency offset due to the Doppler shift over the propagation 
* channel (a fimction of speed, v, with respect to the BSTs and Direction 

Of Travel (DOT), q), with respect to Northing). 

Xj = TOA of a signal pi(t- j ) (a fimctidn of the range Rj), 

Ati is the overall group delay through the i*^ BST to the MR antenna 

(usually a fimction of frequency), 
Wi(t) represents the Additive White Gaussian Noise (AWGN, complex) due 
to thermal noise at the MR (a fimction of temperature and bandwidth), 
MPi(t) represents all multipath components (complex) at MR antenna (a 
fimction of the environment and of the elevation of the antennas), and 
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^(1) represents low-pass equivalent interference (both in-band and out-of- 
band) at the MR. 

Total RF Reception: In IS-95, a number of pilot signals, sync signals, paging 
signals, and forward traffic signals, are transmitted over the same frequency 
band in the forward link. In GPS, a number of signals are transmitted from the 
satellites over the same frequency band. In other words, in a CDMA system 
^(t) represents aU CDMA signals excluding q(t). 

Baseband Reception: Several IF stages down-convert the received RF signal. 
rj(t), to an analog baseband signal, fjCt), which is equivalent to performing the 
following operation: 

f,(t) = Koic^-^ m e-^^"^'-*^>' e*- }= Re{^(t) e^^"^- ' e'^ (8) 

whore 

• is the frequency error between the Local Oscillators (LOs) in the 
MR receiver and the carrier frequenQr f^; and 

• ^ is the carrier phase in Ae Local Oscillators (LOs) of the MR. 

Ft-om equations (5), (6) and (8), one can refer fe>; 

• • Vi = 0i - Co + Q" as the Phase Of Arrival (POA) of the received signal, 
fj (t)., from the i* BST at the MR; and 

• % = 5f; - A^, + Af;" as the Frequency Of Arriyal (FOA) of the received 
signal, ii (t) , from the i'' BST at the MR. 

ni pitfll Reception: Finally, an Analog-to-Digital (A-to-D) Converter converts 
the baseband sigioal, f(t), to a digital signal, r(m), which is equivalent to 
performing the following operation: 

?(ni) = f(t)U^T+q(m) (9) 

where 

• AT is the sampling interval in the MR; 

• m is an integer; and 



SUBSTITUTE SHEET (RULE 26) 



wo 01/65271 




PCT/CAOO/00224 



30 

• q(m) is the quantization noise which depends on the number of bits in 
the A-to-D. 

Sources of Location Information in equ ation (S): 

There are several sources of information where one can extract an independent 
equation from equation (5) regarding the static location of the MR, i.e. 
regarding the coordinates, (x,y), of the MR: 

1. Rp. 

2. \j/i,and 

3. Pi(t-Ti+X,.i-Ati). 

In addition, there is one source of information where one cani extract an 
independent equation from equation (5) regarding the kinematic location of the 
MR, i.e. regarding the velodty, v , namely: 

Starting with &e static information, we have: 

jL_Rr" : In the static case, tiie first source, RJ"*, can be used to estimate Ri- 
However, Rp is umcliable due to jthe nature of the radio diannel whCTfe RF 
shadowing and flat feding can deviate the value of "n" in "Rr"" from 2 (for 
Line Of Sight (LOS)) to 4 (for an urban environment, as shown in Hata, M., 
"Empirical Formula for Raio Propagation Loss in Land Mobile Radio 
Services," IEEE Transactions on Vehicular Technology. Vol. VT-29, No. 3, 
August 1980), or even 6 in heavy urban environments. Hence, unless LQS is 
guaranteed between the MR and the i* BST, R^" can only offer an accuracy of 
a few kilometers for Rj. In the case whai LOS is guaranteed (e.g. in a flat 
eanvironment such as on the sea or over a lake), G; has to be known to the i 
BST and 'V can be chosen to be 2. 

2. The second source, the POA 

H/i = { -2n{ 5^ +Af; )(Ti - T„_i + Ati ) + ACi }mod 2n (10) 
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can be used to estimate the TOA, X;, where A^i = Ci - Co- Tl^e POA, v|/;, is related 
to the range, R^, between the i* BST, and the MR through equations (10) and 
(11): 

Ri = (Xi-x„.i)c (11) 
where c is the speed of propagation. 

Since the time of transmission, x^j, in (11) is unknown to the MR, one can 
cithen 

L estimate it, then use the estimated value in TOA positioning, or 
2. remove it using TDOA positioning. 
Both TOA and TDOA positioning are explained later in the patent 

In TOA positioning, a possible method to estimate is to estimate the Round 
trip Delay (RTD) between: 

• the transmission from the MR to an active BST and 

• the reception by the MR of the response of the active BST. 
In Has case, 

2Ri = (2Ti.XMR)P (12) 

where Tmr is the time of transmission from the MR to the B.St and 2Ti is the 
TOA of the response of the BST to the MR. A drawback in estimating the RTD 
from equation (li2) is that equation (12) assumes that the BST is able to 
trai&mit a response back to the MR. This assumption makes the handset-based 
WLS depend on the functionality of tiie BSTs. This implies that the handset- 
based WLS is network dependent Nonetheless, we will carry on using such a 
conc^t in order to be able to aid a handset-based solution with a network- 
based solution. 

Furthermore, in TOA positioning, a drawback in using in equation (10) to 
estimate is the existence of a "modulo 2n"" operation in (10) which impUes 
that there can be a large integer ambiguity in the TOA estimation. 

In TDOA positioning, a pair of TO As, and Xj, is needed per independent 
equation. In this case, .the amount of ambiguity is a function of the baseline 
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between the two transmitting antennas, i.e. between the i*** BST and the j*** BST 
with respect to the wavelength corresponding to the frequency: + 5£j H-Afj in 
(10), as shown in equation (13): 

- T,- = (To r -Coj) - {(Vi - v|/j) - (ACi -Ag + 271 n}/27cf, - (At, - At,.) + extra term 
(13) 
where 

• Xj corresponds to the TOAofri(t) at the MR; 

• Tj. corresponds to the TOA of rj(t) at the MR; 

• i corresponds to the TOT of pi(t- x^^; 

• T^j corresponds to the TOT of Pi(t- x^,j); 

• n is the integer ambiguity; and 

• the extra term in (1 3) is due to the existence of 

• Af; m (10) which depends on the error in the LOs for the i*** 
BS relative to f;.; and to 

• 5f; in (10) which dqpends on the velocity v of the MR and 
the wavelength X corresponding to ^. 

The term (x^^ - x^^).is known to the MR and can be included in equation (13). 

In the iS-95 and J-STP-008 standards, the GT frequency, + A^, can deviate 
by up to O.OSpart per million (i.e, by lip to 95Hz for a l900MHz carrier 
frequency). Furthenriore, in cellukr telephony in North America, the carrier 
frequency, f^ can take values around SOOMHz <with a wavelength X of around 
37.5cm) or values aroimd 1.9GHz (with a wavelength X of aroxmd 15.8cm). 
This implies that for a velocity of lOOkm/hr 6^ corresponds to: 

-75Hz > 5fi >75Hz for ^ = SOQMhz, or to 

-150Hz > 5^ >150Hz for 4= 1.9GHz. 

In conclusion: 

1. » 5f; and f^. » Af^, i.e. the extra term in (13) is negligible; and 

2. the amount of ambiguity in (13) can be potentially large unless the 
baseline between the two receiving antennas is small with respect to 
the wavelength, X, corresponding to the frequency: f... We refer to such 
a solution as Solution I. 
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Solution I: When the baseline, d^j, between the antenna at the i'*' BST and the 
antenna at the j*** BST is of the same order of magnitude as the wavelength, X, it 
is possible to generate a solution using either TDOA or PDOA. The TDOA 
solution between the two antennas is a hyperbola while the PDOA solution 
between the two antennas is a line. Both solutions approximate the exact 
solution and coincide asymptotically, i.e. the PDOA line and the TDOA 
hyperbola coincide at infinity. 

Remark: F rom a practical point of view, the baseline, djj, can be of the same 
order of magnitude as X only if the two antennas belong to different sectors of 
the same cellular base station. 

We assurhe that the MR is far from both antennas with respect to d^j. Such an 
assumption implies that: 

1. the transmitted wavefront is planar, i.e. Yj^Yj^Yij, where: 

• Yi is the Angle Of Arrival (AO A) of ri(t) from the antenna of the i*^ BST 
to the MR, 

• Yj is the AOA of r/t) from the antenna of the j* BST to the MR, and 

• Yi J is defined as ttie angle formed between: 

• the line joining the antenna of the i*^ BST to the anteima of the 

BST (which we tefef to as liney); and, 

• the line joining.the MR to the middle of liney between the two 
antCTnas (which we refer to as line^jy) 

in a clockwise manner fix)ni line^^j to liney. 

2, the following equation relates the Phase Difference of Arrival (PDOA), 
(Vi-VjX to the angle Yij: 

— -jsi cos(y^j) = (Vi - )mod 27t (14) 

where: 

• is the Phase Of Arrival (POA) of r|(t) at the MR; and 

• v}/j is the POA of tjCt) at the MR. 

The solution forYg in (14) is 
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(15) 



Yjj =±cos-']-^(M/i-Vj + 27ik) 



where k in (15) is an integer that has to satisfy the following condition: 

_ M(M/.-M/j)niod27r^ 1 (16) 

For example, when djj = X,, we have the following possible solutions for k: 

when (-^i-^Mo^^^^.i^ wehave k = 0.1or2; 

2% 

when wehave k = 0..1or-2; 

2it 

when ^^'"^^^"^"^^" = -1/2, wehave k = Oorl; (17) 
2n 

when ^^^-^^^^^^^"=4-1/2, we have k = Oor-l; 
2n 

^hen ^M^i)^?5^:,o, wehave k = 0.1or-l; 

2n 

Some of the solutions in (17) are trivial. As d^^ becomes larger than X, we h^ve 
more nontrivial solutions. Even though having more than one solution unplies 
ambiguity, it is possible to resolve the ambiguity using Least-Squares 
Techniques as explained later. 

In processing to estimate either or Yi, we have the following sources of 
errors: 

• The effect of noise. w/t>. on y. : 

When using a linear array of antennas to estimate the AOA of an RF 
signal approximated as a planar wavefront (i.e. assuming a distant RF 
soiu-ce), the Cramer-Rao Lower Bound on . the variance, var(Yi), of the 
estimated AOA, yj, is equal to 
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\2c' 



vai 



• iTonc 

where 

• var(yi)[^ ^ is the variance of the estimate y-, of the bearing 
Yi. 

• SNR4^„, = A?/(2a^). 

• Ai is the amplitude of the tone, 

• is the noise variance, 

• M is the number of elements in the antenna array, 

• d is the distance between antenna elements, 

• ^ is the frequency of the tone, and 

• c is the speed of light 

Equation (18) is based on the transmission of a narrowband tone from 
each element of the array at the same phase. When the transmitted 
signal Si(t) is either a wideband signal (e.g. CDMA) or Si(t) has a 
unique frequency (e.g. FDMA), (18) no longer applies directly, but 
can be made to apply after some niathematical adjustments to the 
received signals usmg the Reference Information estimated by flie RR. 

• The effect of noise, w,(tV on Mf^: The Cramer-Rao Lower Bound 
provides a lower bound on the variance, vai^vj)!^^, of the estunate, 
4r . , of the phase vj/j for a tone of frequency due to AWGN : 



(19) 



SNRjL xN5x4jt^f„^ 

Tone 'iTone " 



where 



• varfR- )L is the variance of the estimate, R| of the range R^, 

' •hrone 

• X is the wavelength of the tone, 

• 4, is the frequency of the tone corresponding to X(i.e. 

• SNRi[j.^^^ = A?/(2a^) where SNRj is the Signal-to-Noise 

Ratio at the MR from the i*^ BST. 

• Aj is the amplitude of flie tone at the MR from the i* BST, 

• is the noise variance, 
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• Nj is the number of samples (which is directly related to the 
observation interval Tj through Tf^W^s where ^ is the 
sampling frequency), 

• c is the speed of propagation. 

Once again equation (19) is based on the transmission of a narrowband 
tone from each element of the array at the same phase. When the 
transmitted signal Si(t) is either a wideband signal (e.g. CDMA) or Si(t) 
has a unique frequency 4(c.g. FDMA), (19) no longer applies directly, 
but can be made to apply after some mathematical adjustments to the 
received signals using the Reference Information estimated by the RR. 

• Phase offsets: The phase of the tone 4 is shifted by AQ. When AQ is a 
function of the frequency 4,» it tas to be estimated otherwise the phase 
\|/|is distorted. 

• time delays. At; The propagation delay Xj is affected by At^ which 
represents the systeirn delay tiirough the antenna, cables, filters, 
amplifiers, etc. Once again if At^ is a function of the frequency f^^ it 
has to be estimated otherwise the propagation delay can vary 
significantly from one received signal to another. 

• multi0ath. MP.m: In TDOA, (he mtdtipath MPi(t) is equivaleiit to 
extra delay over the propagation chaimel and has to be either estimated 
and reindyed, or Mtigatted, .W the ease where the phase estimate is used 
to solve for the AOA of the ieceived signal ri(t), the effect of multipatti 
IS to shift the AOA ofx^t} by some amount depending on the AOA of 
MPi(t) and its inagnitude, 

• interference, ^(t^; De;^ehding on the level of interference, X(Ocan 

have a drastic effect on the accuracy of the estimated phase and may 
saturate the RF front end of the receiver. Its effects can be mitigated 
with both analog and digital hardware, as well as adequate software, 
whether it is in-band or out-of-band. 

The practical description describes methods and apparatus to estimate and 
to mitigate its sources of errors. 

3, Pi^t - t, 4- x^ , - At|): The last (and most common) source of information one 
can use to extract an independent equation for the static location of the MR 
is 



SUBSTITUTE SHEET (RULE 26) 



wo 01/65271 




PCT/CAOO/00224 



37 

p,(t-ti + T,,.AO (20) 
which can be used to estimate either or 

• When Xj is estimated directty we refer to the solution as Solution II. In 
this case, either TOA positioning or TDOA positioning can be used. 

• When Yi is estimated from Xi, we refer to the solution as Solution m. In 
this case, AOA positioning is used. 

• When T| is estimated indirectly based on the phases of tones in 
p,(t - Xi + x^ - Ati), we refer to the solution as Solution IV, Once again, 
either TOA positioning or TDOA positioning can be used as discussed 
above. 

Solution n: The most common method for estimating Xi is to cross-correlate 
p.(t . xj + T^^j - Ati) with Pi(t), i.e. to cross-correlate r(t) with Pi(t), In general 
several peaks result from such a correlation particularly when receiving pilot 
signals in CDMA where all transmitted pilots are identical except for a 
relative delay between pilots. In fliis case, the peak corresponding to Tj in Pi(t 
- Xf + Xo i - At^ is. located within a reasonable propagation delay, Xj - % with 
respect to a clock derived from the signal; and can be distinguished from 
other pilots. 

• In TOA positioning, the time of transmission, x^^i , from tiie antenna of 
the i*^ BST to flie MR has to be estimated in order to estiknate Xj - x„ j . 
A possible mefliod is to estimate the RTD between the MR 
transmission and the reception at the MR of the BST response to the 
MR transmission as shown in equation (12). 

• In TDOA ppsitiohingi the time of transinission, x^ j , is estimated using 
ah extra independent equation. This is accomplished by estimating Xj. 
Once again, the most common method for estimating Xj is to cross- 
correlate pj(t - Xj + • Atj) with pj(t), i.e. to cross-correlate rj(t) with 
Pj(t). In CDMA, the resulting correlation peak corresponding to Xj in 
Pj(t - "Cj + 'Coj • Atj) is located within a reasonable propagation delay, Xj - 
x^j, with respect to a clock derived from the signal. The TDOA, Xj - Xj, 
can be used to estimate the relative delay, (xj - x^ 0 - (xj - x^j) since the 
relative TOT, x^y x^j, between the i*^ pilot and the j"^ pilot is known at 
die MR. 
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Solution ni: In order to estimate Yjj from the TDOA solution, Tj - x^, we 
assume that the MR is far from both the antenna at the i^ BST and the 
antenna at the j"^ BST with respect to their basehne, d^j. As mentioned 
previously, such an assumption practically implies that: 

1. Both antennas exist at the same BST. 

2. The received wavefront is planar, i.e. = Yj ^ Yi j, where 

• Yi is the Angle Of Arrival (AOA) of ri(t) from the antenna of the i*^ 
BST to the MR, 

• Yj is the AOA of rj(t) from the antenna of the j*^ BST to the MR, and 

• Yij is defined as the angle foraicd between: 

• the line joining the antenna of the i* BST to the antenna of 
tiie BST (which we refer to as lincy); and, 

• the line joining the MR to the middle of liney between the 
two antennas (which we refer to as linCpjy) 

ill a clockwise maimer from line^^ to liney . 

3. Equation (21) relates the angle y,j to the Time Difference of Arrival 
(TDOA), Ti - tj, as follows 



d. . 

^COS(Yij) =Tj.T3 
O 



(21) 



where 



• dij is tiie distance between the antenna at the i"^ BST and the antenna at 
thej*BST;and 

• c is the speed of propagation. 

The advantage in estimating the AOA, Yij > ^i^g Pi(t - Xj + x^^ - Atj) in 
equation (20) firom equation (22): 



Yij=±cos ^j^(^i-'^j) 



(22) 



instead of using Vj/j in (10) firom equation (23): 
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^-1 



y-. =±cos-'^— ^(M/i-V|/+27ck) 
l2Kd,. 



(23) 



is the fact that there are no ambiguities in equation (22) except for the ± in 
the solution of Yj, while equation (23) can have a number of ambiguities 
depending on the value of the integer k. 

The disadvantages in estimating the AO A, Yj j , using Pi(t - tj + - Atj) in 

equation (20) instead of using xj/; in (10) are: 

L a small baseline is required between the two antennas; and 

2- the noise affects the variance, var(Xj)|p^^j, of T; substantially more 

than it affects the variance, var(Hi^ j ) , of xj/j . 
The second disadvantage can be explained by comparing the Cramer-Rao 
Lower Bound for \{f j : 



with the GRLB for % : 



_^ (24) 



Tone SNR4^xN,x4kV 



c- vartx; )! ^ : = =■ (25) 



s^var(Xi)| 



(«) SNRJ ,.,x4re*BW* 

Assuming that SNRj]^^ xN. in (24) is equal to SNRij^^^j in (25), then the 

difference between (24) and (25) is in (24) compared to BW^ in (25). In 
AMPs the BW is approximately equal to 30KHz while ^ can be equal to the 
carrier frequency which is equal to 800MHz, The ratio between the two 
values in dB is 44,26. In IS-95 and J-STD-008 the BW is approximately 
equal to 1.23MHz while can be equal to the carrier frequency f^ i.e. 
800MHz. The ratio between the two values in dB is 28.13dB. In both cases 
the difference is large. Nonetheless, pj(t - - At^) in (20) is suitable for 
TDOA positioning while xj/j in (10) is suitable for AOA positioning as 
explained here: 

• Using in (10) to estimate the AOA, Yij» we have a lower-bound of 
5x10'^ radians^ for a SNR, SNRL x n., of 30dB and for a distance, d, 

MTonc ' 
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of one wavelength, i.e. 0.375m. Tliis is equivalent to having a standard 
deviation for y, j of 7.07x10"^ radians, which corresponds to a range 
error of 7.11 meters for every kilometer range between the MR and the 
antennas. 

• On the other hand, Pi(t - x-, + x^i - At^) in (20) is suitable for TDOA 
positioning since in this case the baseline between antennas does not have 
to be small. For example, assuming the SNR, SNRil^^^^^ , in equation (25) to 

be 30dB for AMPs, the lower-boimd on the variance, 
var(R|)| =0^ var(Xi)l , of the range R, is 2553meters^ or cquivalently 

R. has a standard deviation of 50.53meters. On the other hand, assuming 
the SNR, SNR l , in equation (25) to be lOdB for IS-95, the range is 

150,69meters^ or cquivalently Rj has a standard deviation of 12.28meters. 

Solution IVt When tj is estimated indirectly based on the phases of tones in 
Pi(t - Xi + - Atj), an integer ambiguity can still exist in the solution unless 
the wavelength of the selected tones is large with respect to the range 
between the MR and the antenna at the i^ BST. The selected tones are 
discussed below where a solution for this integer ambiguity is presented. 

In processing Pj(t - + x„j - Atj) to estimate Xj , we have flic following sources 
of errors: 

• noise. w,(t): The Cramer-Rao Lower Bound provides a lower boxmd on 
the effect of AWGN on the variance, var(Xi)|p^^^j , of the estimate of the 

delay x^: 

var(R:)| =c^var(x:)I ^ r- (26) 

*>Pi(0 

where 

• varfR. )| is the variance of the estimate, Rj , of the range R^, 

* »P.(0 

• BW is the RMS bandwidth of Pi(t), 

• SNR,|^^,^ = 2E,/N„ 

• Ej is the energy of r;(t), 

• No is the noise Power Spectral Density, and 

• c is the speed of light. 
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• The effect of Bandwidth. BW: 
Equation (26) demonstrates that the bandwidth of the signal plays an 
important role in the accuracy of the wireless location system. In IS-95, 
the radio frequency (RF) channels have a bandwidth of approximately 
1.25 MHz which is comparable to that of systems designed primarily for 
location such as Global Positioning Systems (GPS) with a BW of IMHz 
over Standard Positioning Services (SFS) channels and ISM-based 
location systems with a typical bandwidth of lOMHz. In GPS, the location 
system may initially use a conventional sliding conrelator to obtain a set of 
pseudo-ranges (one pseudo-range per satellite) (see e.g. Spilker, JJ., "GPS 
Signal Structure and Performance Characteristics/' Global Positioning 

. System. Volume I The Institute of Navigation, Washington D.C., 1980). 
The pseudo-ranges are then used in midti-lateration to obtain a position fix 
of the GPS receiver. A typical accuracy for a commercial one point (i.e. no 
differential reception) GPS receiver with SPS is aroxmd 30ni RMS without 
Selective Availability (SA). In <iirect proportions, an IS-95 land-based 
location system which initially uses a conventional sliding correlator to 
obtain a TOA estimate of the transmitted radio signal followed, by a 
hyperbolic (differential) multi-lateratibn of all the TOA estimates should 
offer a location accuracy better than 30m RMS asstmiing ho multipath. 
The reason is that the. cpirelatipn J^ctiqn^^^.fi^^ a TO A may be 

estiihateid has a reisolutioii which is'liimted to that pf the Fourier transform. 
The traditional resolution bound on Fourier-based methods is the Rayleigh 
resolution criterion as shown in Haykin, S., Adaptive Filter TTieory/' 2nd 
Edition, Prentice Hall, Englewood Clifife, NJ, 1991. Further processing 
using SR and inverse SR ^gorithms often yields a result with higher 
resplution as shown by Dumont, L.R., et al., "Super-resolution of 
Multipath Channels in a Spread Spectrum Location System," lEE 
Electronic Letters, Vol. 30, No. 19, pp. 1583-1584, September 15, 1994 
and as shown by Fattouche et al., U.S. Patent No. 5,570,305 issued Oct, 
1996, and as shown by Ziskind, L et al., "Maximum Likelihood 
Localization of Multiple Sources by Alternating Projection," IEEE 
Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP-36, 
No. 10, October, 1988. It is proposed to use SR and inverse SR 
algorithms to improve the resolution of the WLS beyond the Rayleigh 
resolution, as opposed to 

Dent, U.S. Patent No. 5,404,376 issued April 4, 1995; 
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Dunn et al., U.S. Patent No. 5,600,706 issued Feb. 4, 1997; 

Khan ct al., U.S. Patent No. 5,646,632 issued July 8, 1997; Ghosh et 

al., US, Patent No. 5,646,632 issued April 16, 1996; 

Krasner, U.S. Patent No. 5,663,734 issued Sept. 2, 1997; and 
Krasner, U.S. Patent No. 5,781,156 issued July , 14, 1998; 
which do not further process the correlation peak. 

• The Effect of the Observation IntervaU T^: 

The Observation interval, Tj , is directly related to the energy in the 
received signal ri(t) at the MR since 

E,=^/h(t)fdt (27) 

0 

It is also related to the steady state response of a filter in such a way that 
the response of the filter reaches its steady state as long as 

T>-^ (28) 

where BW| is the bandwidth of the filter that is processing ri(t). 

• time delays. At: The propagation delay iis affected by Atj which 
represents the system delay tlirougli aiiteima, cables, filters, amplifiers, 
etc. When At; changes from one EST to anptticr independently fi-om the 
RR, it has to be estimated otherwise the propagation delay can vary 
significantly from one received signal to another. 

• interference, ; Depending on the level of interference, Ii(t) can 

have a drastic effect on the accuracy of the estimated phase and may 
saturate the RF front end of the receiver. Its effects can be mitigated 
with both analog and digital hardware, as well as adequate software, 
whether it is in-band or out-of-band, 

• multipath. MP,(tl: The multipath MPi(t) is equivalent to an extra delay 
over the propagation charmel and has to be either estimated and 
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removed, or mitigated. In the case the estimated time delay, Xj, is used 
to solve for the TDOA between the received signal ri(t) and the 
received signal rj(t), the difference in multipath, MPi(t) - MPj(t) , is the 
distorting factor which has to be mitigated. 

• The effect ofMultipath, MP/tt: 
When multipath is considered, the accuracy of the AMPs land-based 
WLS could potentially degrade even further. The cellular frequency 
band is between 824 and 894 MHz and the propagation characteristics 
at these UHF frequencies will have a significant impact on positioning 
by multi-lateration as shown in Parsons D., ''The Mobile Radio 
Propagation Channel, " John Wiley & Sons, New York, 1992. That the 
ranges measured correspond to Line Of Sight (LOS) distances is a 
major assumption made when estimating position by multilateration. 
Although tiie dominant transmission mode in this band is LOS, 
reflectioiis from natural aiid man-made objects as well as diffraction 
arotmd said objects are also possibilities. Multipath and diffraction 
allow die cellular signal to propagate in heavily built up areas as well 
as indoors. However, they also cause the measured ranges to be longer 
than the true LOS distance which introduces error into the 
multilateration process. In addition, the propagation distence at UHF is 
relatively sliort This, allows frequency iieuse in die cellular system but 
limits the number of observables in die multilateration process. For 
instance, in a dense urban environment witii a delay spread of 3 
microseconds (as shown in Hata, M., '^Empirical Formula for Radio 
Propagation Loss in Land Mobile Radio Services," IEEE Transactions 
on Vehicular Technology, Vol. VT-29, No. 3, August 1980) multipath 
causes the location accuracy to degrade to more than 1400m RMS. 
Once again, the reason for this is that the correlation ftmction from 
which the muitipath may be estimated has a resolution which is limited 
to that of the Fourier transform which implies that any multipath within 
such a resolution is unresolvablc using traditional methods. Further 
processing using an inverse SR algorithm often yields a result with 
higher multipath resolution as shown by Dumont, L.R., et al., " Super- 
resolution of Multipath Channels in a Spread Spectrum Location 
System," lEE Electronic Letters. Vol. 30, No. 19, pp. 1583-1584, 
September 15, 1994 and as shown by Fattouche et al., U.S. Patent No. 
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5,570,305 issued Oct., 1996, and as shown by Ziskind, 1. et aL, 
"Maximum Likelihood Localization of Multiple Sources by 
Alternating Projection," IEEE Transactions on Acoustics, Speech, and 
Signal Processing, Vol. ASSP-36, No. 10, October, 1988. Another 
approach to multipath resolution is due to Morley, G.D. et al., 
"Improved Location Estimation with pulse-ranging in presence of 
shadowing and multipath excess-delay effects," Electronics Letters, 
VoL 31, No. 18, pp. 1609-1610, 3lst Aug. , 1995, It is proposed to use 
SR and inverse SR algorithms to better resolve the multipath 
components in the received radio signal, as opposed to 

Dent, U.S. Patent No. 5,404,376 issued April 4, 1995; , 
Dunn et al., U.S. Patent No. 5,600,706 issued Feb. 4, 1997; 
Khan et al., U.S. Patent No. 5,646,632 issued July 8, 1997; 

Ghosh et al., U.S. Patent No. 5,646,632 issued April 16, 1996; 
Krasner, U.S. Patent No. 5,663,734 issued Sept. 2, 1997; and 
Krasner, U.S, Patent No. 5,781,156 issued July , 14, 1998; 
which do not fiitthCT process the correlatioti peak. 

Ghbices of Tones for Solution IV: 

irhere are existing tones and there are generated tones. Both types of tones 
cany the MR Location Mbimation in their phases (mod 27t). 
Existing Tones in thaperiodic pilot sequence: 

When the received signal, ri(t), is periodic (e.g. the pilot signal in IS-95, 
the CA-<x>de in GPS)^ it can be approximated using Fourier series, or 
equivalenfly, a number of tones can be used to approximate the pilot 
sequence: 



where 

• 4i is the fundamental frequency of the pilot sequence; 

• Tj is the propagation delay at the antenna of the i^ BST; and 

• c„ is the complex Fourier Coefficient of the pilot sequence. 
Using equation (30) below one can estimate the TO A, t,-, as follows: 

(t-Xi) = {(phase of n*^ term in (29)) modulo 27t}/(27cnf J (30) 



CO 



Si(t-Xi)= Sc„exp[27mf„(t-Xi)] 



(29) 
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Geperated Tones using a nonlinear operation on the pilot scq uenro* 
When the received signal, ri(t), is cycle-stationary, it generates a tone 
when a nonlinear operation is performed over it. The phase of such a tone 
carries location information regarding the MR. Most digital signals are 
cycle-stationary and the generated tone is usually at the symbol rate 
(except in CDMA where the generated tone is the chipping rate). 

In the case when the BSTs are AMPS/TDMA/GSM-based, taking 
advantage of the (ryrclo-stationarity nature of the received signal is 
particularly important since in this case, it is not important to have a clean 
replica of the received signal. After performing the nonlinear operation on 
the received signal, a tone is generated of known (or estimated) frequency. 
This implies that the duration of observation can be extended as long as 
there is a signal to be received, i.e. as long as a call is in progress 
regarding a CT. This leads to an improved SNR for the received, signal as 
explained in the CKLBs shown above. 

The practical description introduces methods and apparatus to estimate Xj and 
to mitigate its sources of errors. 

4. f;; In the kinematic case where the MR nioves with respect to the BSTs, the 
speed and Direction Of Travel (DOT) of the MR is of interest. The only source 
of infoimation one can use to extract an independent equation for the speed and 
DOToftheMRis. 

^ = 5^..A^, + A£; (31) 

which can be used to estimate the Doppler shift 5f^ . This is explained in Figure 
4 which assumes a three antenna system: the antenna at the m^ BST with 
coordinates (x„,yni,2jj, the antenna at the k**' BST with coordinates (x^.y^jZ^ 
and the antenna of the n*^ BST with coordinates (x„,y„,zj. In this case, the 
FOAs: f^y ffc and 4 are related 

• to the Direction Of Travel (DOT), (p, relative to Northing (clockwise), 

• to the relative speed of travel, v, with respect to the BSTs, and 

• to the frequency offsets A^, as follows: 
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f;„ = f, + 6f^- Af, + Af;,= f, + vA cos((p- AJ - Af, + Af;„ (32a) 
f, = f, + 6fk - Af, + Af, = f, + vlX cos((p- A,) - Af, + Af, (32b) 
f„ = fc + 5£;.- Afo + Af;, = £e + vA, cos((p- A„) - Af, + At;, (32c) 

where 

• is the clockwise angle from Northing to the line formed by (x„,y J and 
(x,y), 

• A,, is the clockwise angle from Northing to the line formed by (x,,,yj and 
(x:,y), 

• A„ is the clockwise angle from Northing to the line formed by (x„,yj and 
(x,y),and 

• (x,y) are die 2-D (horizontal) coordinates of the MR. 

In other words, there are three equations (32a, b and c) with 6 unknowns: <p, v, 
A4, Affc, Af; and A^. This is a problem that can be resolved if the frequency 
ofifsets: A4, Afj,, A4 are estimated by the RR and removed. Then we have two 
equations with two unknowns: v and (p after using Frequency Difference Of 
Arrival (FDOA): 

4 - 4 = vA, cos(q)rAJ| - vA, cos(q)-AO (33a) 
= v/X cos(q)rAfc) - vA. c?os(<p-AJ (33b) 

We refer to such a solution as Soliitipn V. Note that in Solution V the position 
of the MR (x,y) must be known (or estimated) prior to estimating the speed and 
velocity of the MR, iii order to be able to Imow A^, A^ and A„ in equations (33). 

Iii processing f; to estimate 6fj , we have the following soiuces of errors: 
• noise, W|ft): The Craiher-Rao Lower Bound provides a lower bound on the 
variance, var(5fi)L , of the estimate of the frequency 6fl of a tone due to 

■Tone 

AWGN: 



/ar(5fi)L 

IToi 



12 



Tone SNR;!^^^ X Ni X (N;^ - 1) X 471^ 

where 

• var(Sf.)| is the variance of the estimate, 5fj , of the frequency Sf^, 

' Hrone 



(34) 
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• SNR-L = A? /(2(jI) is the SNR of the tone at the MR from the antenna 
of the i* BST, 

• A,- is the Jtmplitude of the tone at the MR from the antenna of the i*^ BST, 

• is the noise variance, and 

• Nj is the number of samples (which is directly related to the observation 
interval Tj). 

• frequency offsets, Af,; As mentioned above, there are three equations (32a- 
32c) with 6 unknowias: <p, v, A£;„, Af^, Af;, and A^. This problem can be 
resolved if the frequency offsets Af) are estimated by the RR and removed. 

• multipath. MP.m: The effect of the multipath MPi(t) in this case is to add 
some Dpppler shift to 6^ due to dynamic (non-stationary) reflectors such as 
cars, buses, trucks, etc. 

• interference. ; In this case, ^(t) plays the same role in distorting the 

estimate of the frequency 5^- as in distorting the estimate of the phase G| 
except fliat its effect can be reduced more significantly. 

• Error in Estimating fx.v>: In equations (33) it is assumed that the position, 
(x,y), of the MR is known prior to estimating its speed and DOT. This is 
usually not true and (x,y) needs to be estimated. fir^t The estimation of (x,y) 
is imperf^^ ioaplying that it will contsun errors that can affect the estimation 
of the speed of tiiie MR^^^^ 

The practical description explains methods and apparatus to estimate 64i» 54 
and 54 and to mitigate sources of errors. 

Positioning Technologies: 

1. TOA/Range Positioning (Circular Multilateration) 

In a 2-D (horizontal) location system it is possible to estimate the position of a 
MR from the range of tiie MR to at least two BSTs, The MR horizontal 
position estimate is simply the intersection of two horizontal circles with radii 
equal to the ranges and centered at the BSTs. This method of positioning may 
therefore be called circular multilateration. Ranges may be calculated by 
subtracting the known time of transmission of signal s(t) from the measured 
TOA, Tj, of signal ri(t). There are three possible solutions: 
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1. When the Time of Transmission, x^, is known, we refer to the positioning 
technique as Range positioning. 

2. When the time of transmission, x^,, is unknown, and to be estimated, we refer 
to the positioning technique as TOA positioning. 

3. When the time of transmission, x^, is unknown, and to be eliminated using 
TDOA, we refer to the positioning technique as TDOA positioning (which 
is discussed in the following two sections). 

In Range positioning, one way of determining the time of transmission, x^, is to 
use the RTD between a MR and a BST with respect to the round trip delay 
from tiie time the MR transmits a signal to the BST to the time it receives the 
response of the BST to the MR transmission. In this case, a minimum of two 
independent equations is required to solve for the two unknowns x and y. If 
two ranges arc estimated from TOA measurements of two different BST's, then 
there are two solutions for MR location at the intersections of two circles 
centered at the two BST*s. This ambiguity can be resolved by including a third 
range measurement from another BST, or by (approximate) angle estimation 
firoin one or more sites by examining signal strengtii or phase difference of 
arrival. 

fri TOA positioning, there are now three unkrro\sa;s: x, y and x^. TOAs from 
three BSs are required and the equation corresponding to the i*^ BST is 

x,-Co--V(x-Xir+(y-yif =0 (35) 

C 

where 

• Xj is the Time Of Arrival of signal Xi(i) at the MR from 
the antenna of the i"^ BST, 

• Xj, is the time of transmission of signal s(t) from the i* 

BST, 

• (x, y ) is the best known 2-D position of , the MR, 

• (xi>yi) is the known 2-D pdsition of the antenna at the 
i^ BST where i=l, .,.,N. 
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When more than the minimum number of BSTs are available, as explained 
above, redundancy is said to exist. If a redundant set of measured TOAs 
contains errors, the TOAs must be adjusted in order to obtain a unique solution 
to the unknowns. The amount by which each TOA is adjusted is called the 
residual. Several methods of adjusting the observed according to some 
optindzation criteria are available including: 

• General Minimum Variance Unbiased Estimation, 

• Best Linesu: Unbiased Estimation 

• Maximum Likelihood Estimation, 

• Least Squares Estimation, 

• Method of Moments, 

• General Bayesian Estimation, 

• Linear Bayesian Estimation, 

• Kalman Filtering, etc. 

It is preferred to adjust the observed TOAs such that the sum of the squares of 
the residuals is a minimum. This is Least Squares Estimation, 

Effect of Geometry. on TOA Positioning 

In a location system, geometry (that is the relative positions of the BSTs with 
respect to each other and the MR to be positioned) plays an important role as 
shown in equation (36) 

std(position) = DOP x std(Ri) (36) 

where 

• std(Ri) is the standard deviation of the range estimate 
R. obtained from equations (19) or (26), and 

• DOP is Dilution Of Precision which is a measure of 
geometry. 

In a 2-D location system, it is common to refer to the Horizontal DOP (HDOP) 

which is defined as 

HDOP = VeDOP^+NDOP^ (37) 

where 

• EDO? is the East DO? and is defined as the square root 
of the element in the row and 1^' column of Ctoa. 
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• NDOP is the North DOP and is defined as the square 
root of the element in the 2"*^ row and column of 

Ctoa> 

• Ctoa ^ unsealed matrix defined as 

• Cj is the unsealed measurement covariance matrix (the 
identity matrix of appropriate dimension), 

• A is the design matrix for the model of equation (35) and 
is defined as 



A = i 

Q 



x-x 



X-X2 



1, _ 



y-Yi 
d, 

y-y2 



-c 



— c 



(38) 



and di is the best deriy€?d distance (range) between tiie MR and the antenna at 
thei'^BST. 



Weighting of the TQA Observations 

In the position estimiatiori process, not all TOA observations need cany the 
same weight. TOA observations thought to be more reliable may be weighted 
more heavily than those deemed less reliable. This is accomplished through 
the observation covariance matrix. The inverse of Cj is the weight matrix. 
Larger values on the diagonal of C,"' correspond to heavier weighting for the 
corresponding TOA observations. In the context of cellulsir telephone 
positioning, RSSI from each BST is one method of assigning weights to the 
TO As. A high RSSI from a BST implies a reliable TOA. This is due to two 
facts: 

1 . The RSSI from a BST usually consists of received signal power + 
received noise power. The noise in the receiver is mainly thermal 
noise which is a function of bandwidth and temperature. When two 
BSs have comparable temperatures and bandwidths, the received noise 
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power is approximately the same in both BSTs. Thus, a high RSSI 
implies a high received signal power, which in turn implies a high 
SNR- This is desirable, 
2. Furthermore, a higher RSSI usually implies less shadowing than a 
lower RSSI, which in turn implies less multipath. This is also 
desirable. 

Blunder Detection in TOA Positioning 

Blunders are gross errors in the TOA observations. They may be caused by 
large signal level fluctuations due to either flat fading or sudden in-band 
interference. If unremoved, blunders cause large errors in the estimated 
position. It is possible to detect observations containing blimders by observing 
the misclosure of each observation during the iterative Least Squares process. 
Misclosure is defined as flie value of the positidn model (35) given the best 
available position estimate (x,y). Observations containing blunders will 
generally have much larger misclosures than observations not containing 
blunders. When an observation is detected as having a blunder it may be 
removed firom the position estimation process. 

We attempt to reduce the efifect of geometry^ Lc. reduce the value of HDOP in 
(37), by allowing the MR to monitor a large number of BSTs. In cellular 
commxmicatiQns, frequency reuse and flat fading are common occurrences. 
Therefore, increasing the nimiber of monitoring BSs (which probably reduces 
HDOP) generally increases std(Ri). As a result blunder detection crucial as 
part of the method of minimizing positional error by maiximizuig the number of 
BSs (and hence reducing HDOP) without iricurring a large penalty on std(Ri). 
In cases where blimders are not detected, poor SNR at some BSs may cause the 
inclusion of measurements from thoise BSs to increase the std(range) more than 
they reduce HDOP- With this trade-off in mind, the number of BSTs used to 
locate the MR is optimized such that std(p6sition) in (36) is minimized. 

Blunders may also be detected by a statistical analysis of the observation 
residuals computed from the Least Squares process as shown by Vanicek, P., 
Krakiwsky, E., ''Geodesy: The Concep/^," North-Holland Publishing Company, 
Amsterdam, 1982. The residual of each TOA observation may be standardized 
by its own estimated standard deviation such that the entire set of residuals is 
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assumed to belong to the normal distribution with zero mean and unit standard 
deviation. If this hypothesis is correct, the standardized residuals should fall 
within some specified confidence region for the standard normal distribution. 
An observation whose residual is flagged is suspected of containing a blunder. 

2, TDOA Positioning fHyperboli o MuUilateration^ 

Instead of estimating the unknown time of transmission, it is possible to 
eliminate it This is accomplished by dififerencing tOAs firom two different 
BSTs. Since the time of transmission fi^om bofli BSTs is known at both the RR 
and MR, it can be eliminated fi^om the resulting TDOA (Time Difference Of 
Arrival). It can be shown that the locus of points for which a particxilar TDOA 
is valid corresponds to a hyperbola. The side of the hyperbola on which the 
MR must lie is known by the sign of the TDOA. Given TOAs firom three BSs, 
two independent TDOAs may be formed The intersection of the two 
corresponding hyperbolas estimates the position of the MR- This mettiod is 
coimnonly referred to as hyperbolic multilateration. 

The 2-D positioning model for hyperbolic multilateration is 

Ar^^4Vtc-x,f + (y-yif-*:^#:^x,y+(y-y,f=.0. (39) 

where 

• Ax{j=i;i-Tj, i.9s^j,and 

• Tj is the Time Of Aiiival pf signal rjCt) at the antenna of 
thei^BST. 

Effect of G^metrv on TDOA Positioninip 

Geometry affects TDOA positioning as well. The HDOP is again calculated 
from (37) where the desigii matrix is now 
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x-x, 
+ 


y-y2 








X-X3 


x-x, 


y-yj 






d3.n 




x-x, 


y-yN 




d, 


dN 



(40) 



d, 



where N in (40) is the number of BSTs. Note that the TOA from the antenna of 
the first BST is subtracted from all other TOAs. 

The imscaled observation covariance matrix is 



2 1 - 1 

1 1 : 

I 1 1 

1 1 2 



(41) 



Weighting of ffieTDQA Observations 

Wei^&ig of the ITDOA observatioris is possible. However^ because the 
TDOAs are a fimction of two TOAs, flie inethod of assigning weights is no 
longer straightforward. 



Blunder Detection in TDOA Positionihg 

Blunder detection may also be performed in TDOA positioning, Howefver, in 
this context, misclosures and residuals correspond to TDOAs. Therefore, a 
failing misclosure or residual may be due to a blunder in either of the TOA 
observations from which flie TD.OA is derived. It is not always possible to 
isolate the offending BST. 

Multiple Solutions in TDOA Positioning 

Two hyperbola halves, formed from two independent TDOAs, may intersect 
twice. This results in two mathematically correct solutions. This is 
particularly true in the case of positioning MRs where short distances and poor 
geometry are commonplace. 
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Solution bifurcation (the existence of two solutions to the positioning 
equations) is most often a concern in the exactly detemnined case. For the 
exactly determined case, the existence of two solutions can be detected with the 
method given in Chaffee, J.W. et al., "Bifiurcation of Pseudorange Equations," 
Proceedings of the 1993 National Technical Meeting, San Francisco, 
California, January 20-22, 1993, The Institute of Navigation, Although 
originally intended for the detection of bifurcation in the GPS (Global 
Positioning System), this method is eiqually applicable to the case of MR 
positioning. 

When more than two TDOAs are available in the 2-D positioning case, the 
probability of exact solution bifurcation is extremely small. It is, however, 
possible for bifurcation to exist for a subset of two TDOAs. In such a case, the 
second solution may affect the final solution obtained when using all available 
TDDAs. 

In any case, when solution bifurcation exists, the iterative Least Squares 
position estimation algorithm may converge to eiflier of the solutions. The 
solution converged to is a function of the initisil starting position u^ to begin 
the iterative Least Squares process. To converge to the solution conre^onding 
to the actual position of the MR, th& initial position used to begin Least 
Squares must be relatively accurate. Given no a priori information about the 
location of the MR, a. closed-form position estimation algorittuh, using the 
TOA or TDOA observations, is the only choice. 

A niraiber of closed-form positioning algorithms have been developed. 
Examples are, spherical interpolation (Smith, J.O., et al., "Closed-Forai Least- 
Squares Source Location Estimation from Range-Difference Measurements," 
IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP- 
35, No. 12, Dec. 1987, pp. 1661-1669), the method of Schau and Robinson 
(Schau, H.C., et aL, "Passive source localization employing intersecting 
spherical surfaces from time-of-arrival differences," IEEE Transactions on 
Acoustics, Speech, and Signal Processing. Vol. ASSP-35, No. 8, Aug. 1987, 
pp. 1223-1225), Bancroft's method (Bancroft, S., "An algebraic solution of the 
GPS equations," IEEE Transactions on Aerospace and Electronic Systems, 
Vol. AES-21, No. 7, Jan, 1985, pp. 56-59), the method of Chan and Ho (Chan, 
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YT., et al., "A Simple and Efficient Estimator for Hyperbolic Location," IEEE 
Transactions on Signal Processing, Vol. 42, No. 8, Aug, 1994, pp. 1905-1915), 
and LOCA (Location On the Conic Axis) by Schmidt, R.O. "A New Approach 
to Geometry of Range Difference Location," IEEE Transactions on Aerospace 
and Electronic Systems. Vol. AES-8, No. 6, Nov. 1972, pp. 821-835. 

LOCA is used in the preferred embodiment of the invention. LOCA is the 
mathematical dual of hyperbolic trilateration. The fimdamental theorem of 
LOCA states that TOA differences from three BSTs of known location yield a 
straight line of position. This straight line is the major axis of a conic. The 
three BSTs lie on the conic and the MR, the location of which is being 
estimated, lies at one of the foci of the conic. In the case of redundancy and 3- 
D positioning, LOCA is expanded into Plane Intersection as shown by 
Schmidt, R.O:, "A New Approach to Geometry of Range Difference Location," 
IEEE Transactions on Aerospace and Electronic Systems, VoL AES-8, No, 6, 
Nov. 1972, pp. 821-835. 

Plane Intersection is equally applicable to the case of redundancy in 2-D 
positioning. 

In LOCA, three conies are possible: an ellipse, hyperbola, or parabola. The 
conic of interest depends on the relative positions of the BSTs and the MR to 
be positioned. Each of these conies has two foci (for the parabola one focus is 
at infinity) and therefore two possible solutions. In the case, of an ellipse, flie 
correct fbciis is determined from the signs of the TDOAs. In the case of the 
parabola, the solution is obvious. When the conic is a hyperbola, however, the 
correct focus cannot be determined from the given TDOAs. Chaffee's method 
of biftircation detection in Chaffee, J.W. et al., "Bifurcation of Pseudorange 
Equations," Proceedings of the 1993 National technical Meeting, San 
Francisco, California, January 20-22, 19.93, The Institute of Navigatioui will 
indicate biftircation for the case of the hyperbola and no biftircation for the 
cases of the ellipse and parabola. With four BSs, a second conic is obtained. 
The intersection of the major axes of the two conies is the MR position 
estimate. Hence, there is no ambiguity and biftircation is, therefore, not of 
concern when using Plane Intersection in the presence of redundancy. 
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If the TO As contain large errors, LOCA or Plane Intersection may give a very 
poor solution. When used as the initial position for Least Squares, this poor 
solution may cause Least Squares to diverge. In that case, the best available 
course of action may be to average the coordinates of the participating BSs and 
use that as the initial position for Least Squares. Another alternative is to use 
the approximate coordinates of the BST with which the MR is conimunicating, 
as the initial position. 

3, Hybrid TDOA Positioning (Circular Multilateration with TDOAsl 
In some cases it is desirable, or even necessary, to use TDOAs as opposed to 
TOAs. However, as explained above, it is not always possible to isolate 
blunders when using TDOAs. It is possible, however, to constmct a 
positioning model that uses TDOAs but gives residuals and misclosures for 
individual TOAs. 

To dp this, one BST is chosen as the reference. The assumption is then made 
that the TOA at the reference site is equal to the time of transmission, x^. 
From this point on, the development follows that of circular multilateration 
except that the TOA of the reference BST is subtracted from all TOAs 
(including that of the reference BST) instead subtracting the time of 
transmission. The positioiung model is then 



At„-FiR„-iV(x-x.)'+(y-yi)'=o (42) 

C C 

where 

• AXj, =ti-Ti, i = l, --,N and 

• Ru is the range from the antenna, of the reference (first) 
BST to the position of the MR. 

The model of (42) has the same structure as that of circular multilateration 
(equation(35)) except that the third unknown is R^ and not the time of 
transmission. 

The model of (42) gives N observations, and therefore N residuals and 
misclosures, for N BSTs. The N observations, however, are a 0 and N-1 
TDOAs. Tlie structure of the equations is that of TOA positioning whereas the 
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observations are that of TDOA positioning. Hence, the name Hybrid TDOA 
Positioning. 

Effect of Geometry on Hybrid TD OA Positioning 
TTie design matrix for hybrid TDOA positioning is 



x-x, 


y-Yi 




d, 


X-Xj 












dN 


dN 



(43) 



Note that (43) is the same as (38) except for the third column which 
corresponds to the third unknown. 

The unweighted and unsealed observation covariance matrix, C,, is the same 
as fliat .for TOA positioning, &e identity matrix of dimension N. 

DOPs m^y be calculated, as before, from.^^C,~*A]' . 
Weiyhtipg of the Hybrid TDOA Oh?5ervations 

Although the observations arc, strictly speaking, TDOAs, they may be treated 
as obsen^tions corresponding to the individual BSTs. That is, the first 
observation, which is always zero, may be viewed as the observation of the 
reference (or first) BST. The second observation, consisting of the TOA of the 
reference BS subtracted from the TOA of the second BST, is treated as the 
observation of the second BST, and so on. Therefore, the N diagonal elements 
of C, may be weighted for the individual BSTs. 

Blunder Detection in Hvbrid TDO A Positioning 

In like manner, N independent misclosures and residuals, one for each BST, are 
available for analysis. Blunder detection may then be performed as in the case 
of TOA positioning and blunders in individual TOA observations may be 
detected. 
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Multiple Solutions in Hvbrid TDQ A Positioning 

Since the observations used are actually TDOAs, this method of positioning is 
subject to solution bifiircation. The method of dealing with multiple solutions 
discussed under hyperbolic multilateration is equally applicable here. 

4. AOA Positioning: 

Given unambiguous AOA measurements, only two unknowns exist for position 
estimation - the 2-D coordinates of the MR to be positioned. Therefore, a 
minimum of two independent AOA measurements are required. Again let 
(x,y) be the 2-D coordinates of the MR and (xj.yj) the 2-D cqordinatcs of the 
antenna of the i* BST. At the antenna of the i* BST, the AOA, A,, is 
measured, where: 

• Aj is the clockwise angle from Northing to the line joining the MR to the 
antenna of the ith BST (similar to Aj in equations (32) and (33); not to be 
confiased with fi in equations (15), (16), (21), (22) and (23)). 

The mathematical.positioning model is&en 

(x-Xi)cosAj -fef-yi)sinAi =0 (44) 

When more than two independraif AOA observationis are available. Least 
Squares may be used to obtain a unique solution. 

fcffect of Geometrv on AOA Positioning 

The relative positions of BSTs with respect to each other and the MR to be 
positioned are also important for AOA positioning. DO? may again be used to 
quantify geometry- From the 2-D positioning model (44), the design matrix is 
found to be 



A = 



cosA| 


-sinAj 




cosAj 


-sinAj 


(45) 




-sin An 
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The unsealed and unweighted observation covariancc matrix, C,, is the identity 
matrix of dimension N. The DOPs can again be found from [a'^C,"* a]" . 

Weighting of the AOA Qbservations 

The individual AOA observations may be weighted through the observation 
covariancc matrix C,. Those AQA observations deemed more reliable are 
assigned a smaller variance, or conversely, a heavier weight. 

Blunder Detection in AOA Positioning 

As in the case of tOA or TDOA observations, AOA observations may contam 
blunders. When redundant observations are available, statistical testing of the 
observation residuals and misclosures may be used to detect observations 
which contain blunders. 

Notes: 

1. As mentioned above, the AOA solution in equation (15) has an ambiguity 
problem which can be resolved eith^ by 

• using TDOA as well as AOA as explained in the next section, or 

• using more than two antennas at the same BS. This is possible in 
sectorized cells since each sector has usually two diversity ant^mas. 
In bflier worSis, in a thre^ sectbr cell a total of six antennas could be 
available. 

2. When the diversity antennas are separated only vertically, the estimated 
AOA is the elevation angle. In this case, it is possible to use both diversity 
antennas since they are. independent from each other, i^e. they offer 
independent observatioiis. Diversity observations miay be combined by well- 
known methods including selection coinbining, maxinial ratio combining, 
co-phasing combiningi equal gain combining, or other methods. 

5, AOA/TDOA Positioning: 

Least Squares allows for the combination of different types of observations. In 
particular, it is possible to estimate 2-D position using both AOA observations 
and TDOA observaitions within a single Least Squares adjustment, 

A combination of the two different observation types is particularly useful in 
situations where only two BSTs are available for estimating the 2-D position of 
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a MR. When two BSTs are used to estimate the incoming signal AOA for 
these BSTs, numerous AOA ambiguities are possible as shown in equation 
(15). This is due both to the nature of AOA estimation with a 1-D linear array, 
and to the fact that the spacing between antennas can be greater than one 
wavelength. Therefore, with two independent AOA observations only, many 
position solutions are possible and there is no way of determining which of the 
many solutions is correct. 

If; however, in addition to the two independent AOA observations, TOA is also 
measiired at each one of two distant BSTs (maybe using Pi(t - + j - Atj) in 
equation (20)), a TDOA can be calculated. This TDOA corresponds to a 
hyperbola side which, in the absence of error, will cross through the 
iaterseption poiiit of two of the many bearing lines derived from the ambiguous 
AOAs. 

When the AOA and TDOA observations include error, no two bearing lines 
and the measured hyperbola side will exactly intersect. For flie measured 
. TDOA and aiiy two particular AOAs, Least Squares will give the position 
solution which minimizes the Sfum of the squares of the residuals. Observation 
residuals aire available since there are two unknowns, flie 2-D coordinates^ with 
three indepie^derifobseryati^^ 1 TDOA)* 

In order to separate the correct AOA pair from the ambiguities, each AOA 
combination is combined with the TDOA observation in Least Squares. That 
combination of AOAs which results in the smallest sum of squares of residuals 
is chosen as correct. The corresponding position solution is used as the MR 
position estimate. 

When combining AOA and TDOA observations in Least Squares, both the 
model for AOA positioning (44) and the model for TDOA positioning (39) are 
used. Without loss of generality, the design matrix A for two independent 
AOA observations and one TDOA observation is 
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A = 



COS A, 



cos A 2 



X — X, 



-sin A, 
-sin A2 



(46) 



where is the best derived distance between the MR and the i*** BST. Note that 
the first two rows correspond to the two AOA measurements whereas the third 
row corresponds to the TDOA measurement The two columns correspond to 
the two imknowns, x and y. Additional AOA and TDOA measurements may 
be included by adding appropriate rows to (46). Note that the 1/c factor seen 
in (40) is missing in the TDOA measurement row of (46). This is done such 
that the imits throughout A are dimensionless. 

The misclosure vector, necessary for the Least Squares adjustment mechanism, 
merely consists of the misclosures of ail observations. The misclosure vector 
corresponding to the design matrix of (46) is 



w = 



(x - Xj )cos A, - (y - y I )sin A, 
(x - X 2 ) cosA2-"(y - y2)sinA2 



c * At,2 - VC'^-Xif +(y^y,)^ -^^li^-^zf HyVif , 



(47) 



where 

• AXl2=Ti-T2. 

Since the AOA and TDOA observations are independent, the unsealed 
observation covariance matrix is the identity matrix of appropriate dimension 
(number of AOA observations plus the number of TDOA observations). 

Effect of Geometry on AOA/TDOA Positioning 

DOP may again be used to quantify geometry. The DOPs (HDOP, EDO?, 

NDOP) can be found from [a^Ci^'a]"^ where A and C, are defined 
immediately above. 
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Weighting of the AOA/TDOA Observations 

The individual AOA and TDOA observations may be weighted through the 
observation covariance matrix C,. Those AOA and/or TDOA observations 
deemed more reliable are assigned a smaller variance, or conversely, a heavier 
weight. 

Blunder Detection in AOA/TDOA Positioning 

Both the TDOA and AOA observations may contain blunders. Statistical 
testing of the observation residuals and misclosures may be used to detect 
observations which contain blimders. 

6. AO A/Range Positioning: 

If the time of transmission or round-trip delay is known, the range, Rj , from 
the MR to the i*^ EST is the observed parameter. In that case, AOA and range 
positioning may also be combined in order to estintiate the position of the MR 
with three BSTs. Without loss of generality, the design matrix A for two AOA 
observations and one range observations is 



A = 



cpsAj 
cosAj 



-sinAj 
-sinAj 

d| J 



(48) 



The misclosure vector corresponding to this design matrix is 



w = 



(x-x,)cosAj -(y-yi)sin Aj 
(x-X2)cosA2 -(y~ y2)sin A2 



(49) 



Since the AOA and range observations are independent, the imscaled 
observation covariance matrix is the identity matrix of appropriate dimension 
(number of AOA observations, plus the number of range observations). 

On the other hand, without loss of generality, the design matrix A for one AOA 
observations and two range observations is 
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A = 



cos A, 


-sin A, 




y-Yi 












d2 



(50) 



The misclosure vector corresponding to this design matrix is 



w = 



(x-x,)cosAi ""(y""yi)sinA, 
fe,- V(x-x.)'4.(y-y.)^ 



(51) 



Effect of Geometry on AOA/Ranye Positioning 

DOP maiy again be used to quantify geometry. The DOPs (HDOP, EDOP, 

NDOP) can be found from where A and C. are defuied 

immediately above. 

Weighting of the AbATRange Observations 

The individual AO A and range observations may be weighted througih the 
obsdrvation covariance matrix C|. Those AOA and/or range observations 
deemed more reliable are assigned a smaller variance, or conversely, a heavier 
weight. 



Blunder Detection in AOA/Range Positioning 

Both the AOA and range observations may contain blimders. Statistical testing 
of the observation residuals and misclosures may be used to detect 
observations which contain blunders. 

7> Speed and Direction of Travel Estimation: 

In a kinematic location system, the 3-D or 2-D location of the moving MR 
must be estimated at various epochs of time. In addition, the doppler shift of 
the signal transmitted from each BST may be estimated as discussed earlier. 
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The equations relating the frequency (including Doppler shift) of the arriving 
signal at the MR from the BST 

1 . to the relative speed of the MR with respect to the BSTs, 

2. to its DOT and 

3. to the frequency offset inherent in each BST 
are given in equations (35). 



The estimation model for speed and DOT is, therefore, 
fj - v/Xcos(9 - Aj) - Afj + Af, = 0 



(52) 



where Afj, the frequency offset for the i*^ BST, is assumed to be known by the 
RR. Given fliat at any particular epoch the position of the MR (x,y) is 
estimated by one of the above methods, the angle Aj in (52) may be calculated 
for each BST. In equation (52) then, the knowns are fj , X, Afj and Aj, (where 
fj is measured or observed), and v, <p and Af„ are the unknowns. Three BSTs 
are required - ^e same number required to estiinate the 2D MR position using 
TO A or TDOA positioning. 

Kffcct of Geometry on Sjieed and Direction of Travel Estimation: • 
Geometry afifects '^e estiipatioh of speed and DOT as well. iPor mstance, it .is 
intiiitively obvious fliat whrant the iS^ is traveling on at line papendicular to the 
line connecting if and a BST, no information regarding the speed of the MR is 
available from observed data at that EST. 

The design matrix for the model of (52) is 



A = 



^^s(2pAj ^3i„(^_^,) _i" 
_cos((t>-Aj Zsin(<p-Aj -I 

_cos(y-AH) l3i„( -1 



(53) 
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The DOPs can again be found from [a'^C,"' a]"* . In this case the DOPs will be 
speed DOP, direction of travel DOP, and frequency offset DOP. C| is the 
unsealed and unweighted identity matrix of dimension N. 

Practical Description of the Invention: 

Forcing a CT to act as an MR: 

In a handset-based WLS, it is desirable to be able to use a cellular BS as a BST 
and to request a CT to act as an MR in order to fulfill the FCC mandate for 
E911 as well as to offer location services to the current 60 million cellular 
subscribers m North America who are expected to grow to 120 million by the 
year 2001. This however might require a change to the CT as well as a change 
in the cellular standards. If the required changes preclude existing CTs, the 
FCC mandate for E91 1 will not be fulfilled tising such a system. Nonetheless, a 
handset-based solution is attractive since it can offer an accuracy comparable 
to the one offered by systems dedicated for location such as GPS and Loran-C. 
This is why it is still acceptable to make changes to the CT regardless of the 
FCC mandate. 

Some chiiges to the CT might require change^ to the standard* We en^hasize 
CT changes while attempting to reduce the changes, to the standard. Three 
levels of CT modification are identified as follows: 

1 . modification 0 which requires no modifications to the CT; 

2. modification I which requires only software modifications to the 
existing CT; and 

3. modification n which requires both software and hardware 
modifications to the existing CT. 

A hardware modification to the existing CT is required when a required change 
in the signal processing algorithms cannot be accomplished in software. This 
may be 

• due to required changes (e.g. a change in tiie RF front end, in the 
antennas or in the mixers) which cannot be performed in software; or 

• due to the fact that the digital hardware is not adequate to perform 
some tasks in software (e.g. when using ASIC chips instead of general 
purpose DSP chips in the existing CTs). 
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It is proposed here to improve modification 0, 1 and II over existing methods. 
In this case, existing methods rely only on pilot signals transmitted from 
cellular sectors to locate the CT. More specifically, existing methods describe a 
typical operation of a CT, which searches for pilots and determines their TO As. 
However, this normal mode of operation of a CT is designed to maximize the 
performance of the CT in receiving data signals, and to maximize the 
communications performance (i.e. capacity of the system, reliability of data 
transmission, etc.). Improvement 0 introduces Location state 0 which the CT 
enters into whenever it is requested to locate itself in modification 0. Location 
state 0 optimizes several design parameters in modification 0 as follows: 
Im provement 0: The most common method to estimate the TOAs, POAs, 
Strength of Arrival (SOAs), and FOAs of received signals at the MR consists 
of correlating a received signal with a clean version of itself- Conventional 
methods cannot resolve corrcilation peaks beyond the conventional time- 
domain Rayleigh resolution. Moreover, conventional methods either choose the 
highest correlation peak above a certain threshold or the peak corresponding to 
the first arrival above a given threshold. Over a multipatfi channel, it is possible 
to assume thai flie first peak above a c^tain threshold corresponds to the 
correct peak, however, aliasing, noise, and interference can cause spurious 
peaks to appear in: time before, the correct one. It is an object of Improvement 
Oto 

1 . optimize the threshold x^sed in selecting the correlation peak; and to 

2. optimize the list of BSs considered for correlation. 

Such an optimizatidn is based on nuiumizing HDOP. The optimized 
parameters: threshold and list of BSs, are then transferred to the CT whenever 
it is requested to enter into Location state 0. In such a state, the CT is forced to 
receive more signals firom more BSTs in order to reduce HDOP. The 
optimization process can be e^lained as follows. la order for ttie CT to receive 
more signals, its threshold has to be lowered. A lower threshold however forces 
the CT to accept weak signals. This produces low SNR signals diat might 
suffer firom multipath and interference. In order to minimize HDOP while at 
the same time reducing the effect of noise and multipath, the optimization 
process consists of solving for flie position of the MR using Least Squares 
weighted by the individual RSSI fi-om each received signal. 
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In modification I, existing methods require a great deal of standard change. 
Instead, Improvement I does not require a standard change. It forces the CT to 
enter into Lrocation state I: 

Tmprovement I: It consists of using the SR and Inverse SR algorithms over the 
correlation function (that is obtained by correlating the received signals at both 
the MR and RR with a clean version of themselves) in order to: 

1. resolve the correlation peaks beyond the conventional time-domain 
Rayleigh resolution; and to 

2. avoid using a threshold in selecting the correlation peak. 

An important feature in the Inverse SR algorithm is the fact that it generally 
produces a number of correlation peaks: 

• a peak corresponding to the direct path signal as well as, 

• a nimiber of spurious peaks caused by noise, aliasing, multipath and 
interference, that can surround the correct peak. 

Instead of using a threshold. Location state I consists of selecting the 
correlation peak that produces the smallest residual when performing a 
weighted Least Squares operation on all possible correlation peaks resulting 
from super-resolving the correlation functions produced from the received 
signals using the inverse SR algorithm. It is also an object of Location state I 
to optimize the list of BSs considered for corpelation based on minimizin g 
HDOP. 

In modification n, existing methods require a great deal of standard change as 
well as Ct change. More specifically, it consists of receiving GPS signals 
partly at the CT and partly at the RR. Instead, several improvements over the 
existing methods are offered in modification BE: 

• It allows the CT to choose between receiving cellular signals, or 
receiving non-cellular signals (such as GPS, Loran-C, etc.), or both; 

• It allows the CT to xise Improvement 0 and/or Improvenient t i.e. to 
enter into Location state 0 and/or Location state 1; 

• It allows the CT to use the carrier phase of the received signals for 
TDOA location; 

• It allows the CT to use AOA for GPS/Loran-C; 

• It allows the CT to locate itself using only two independent BSs, with a 
minimum of two sectors per BS, when using AOA estimation; 

• It allows the CT to use dual-frequency receivers for GPS; 

• It allows the CT to use Pseudolite in GPS; 
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• It allows the CT to use cellular augmented GPS; 

• It allows the CT to use short messaging cither for transmitting its 
Location Information or for receiving the Reference Information, so as 
not to be network-dependent; 

• It allows the CT to combine location information from both GPS and 
cellular signals. GPS alone requires four satellites in order to solve for 
(x,y,z,) and the CT's clock error. Cellular alone require three BSs in 
order to solve for (x,y) and the CT*s clock error. When combining GPS 
equations with cellular equations, we require four equations in total since 
we still have four unknowns: x, y, z, and the CT*s clock error. A possible 
alternative is to replace z by the altitude of the BSs. However, this can 
introduce an error in the estimation for x and y; 

Description of an CT receiving Pilot Signals : 

Here it is proposed to xise an IS-95 receiver which measures the TOA of 
multiple pilot signals for use in estimating the location of said receiven This 
receiver may be obtained through modifications to flie operation of a wide 
variety of CDMA receivers, including CDMA CTs. 

An IS-95 CT is directed to enter location state 0 or Location State I by an 
ord^ received from; 

1 . a remote location, 

2. an ordCT receive from the user of flie CT, or 

3. by otiier predetermined circumstances, such as initialization on power- 
up, or origination of a CT, 

Upon entering location state 0 or I, 

1. the CT assembles a list of pilots to search for (1011) as shown in 
Figure 10. This list may be formed from a combinatiori of the CT's 
Active set, neighbour set, remaining set, and a list of pilots specified 
in a message received by the CT from a remote source. 

2. The process of searching for pilots, and assigning a demodulator 
finger to obtain a refined TOA estimate, is executed for each pilot in 
the pilot search list, hi blocks (1012), (1013), (1014), (1015), (1016), 
(1017), and (1018), P represents one of the pilots from the list. Block 
(1012) initializes P to the first pilot in the search list. 
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3. A searcher element is assigned to search for pilot P in the received 
signal (1013). A window of PN offsets (start and end, or center and 
width) is specified to the searcher. The searcher then computes 
received energy of pilot P at multiple delay offsets within the search 
window. In the preferred embodiment, these delay offsets are spaced 
one half chip apart (0.4069 microseconds). For each delay offset, the 
searcher accumulates received pilot energy over a specified duration 
of time, (preferred embodiment uses 0.5 ms) 

4. The searcher results are examined for evidence of reception of a pilot 
signal with sufficient quality (1014). In the preferred embodiment, the 
ratio of the average pilot chip energy to the total received power 
(called E^) is computed at each delay offset Figure 1 1 illustrates a 
possible set of values for searcher results. In this example, significant 
pilot energy appears near delay offset 7. The pilot EJl^ at each delay 
offset is compared to a threshold value (-15 dB in the preferred 
embodiment)> If the BJI^ at one or more delay offsets exceed the 
threshold, then the earliest such delay ofifset is used as an initial 
estimate of the time of aixival of the pilot If ttie threshold is not 
exceeded at any of the delay offsets, liita it is determined that the pilot 
was not detected. 

5. If pilot P tyas detected, then a demodulator finger is assigned to pilot P 
at the initial estimate of flie tiine of airiyal of said pilot (1015). The 
finger is allowed to track pilot P for a period of time (preferred 
embodiment uses 2 ms), after which 

6. the delay offset of the fingbr is. retrieved and stored ?is the time of 
arrival of pilot P (1016), 

7. After the searcher is finished searching for pilot P, P is assigned to the 
next pilot in the list (1017), (1018). 

8. If the list is exhausted, then the search results are assembled (1019) 
for subsequent processing or transmission to a remote site^ 

Figure 10 illustrates steps performed sequentially. Some steps could be 
performed in parallel. Specifically, the searcher can begin to search for another 
pilot while a demodulator finger tracks a previously found pilot. 
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If the CT's time reference is not sufficiently stable, then one demodulator 
finger can be assigned to one pilot while other pilots are searched and tracked 
by other fingers, in order to maintain a known time reference. 

A demodulator finger may track a pilot using a delay-locked loop, as described 
in U.S. Patent 5,764,687. This results in the finger moving to the peak of the 
arrival of the pilot. Improved performance may be obtained by having a finger 
track towards a predetermined BJl^ threshold on the rising edge of the pilot 
energy versus delay offset. 

Signal processing ; 

The location operation is performed at the Host which can be either located 
within the MR (Option I) or remotely from the MR (Options n and m) (e.g. 
within the RR or a TR). In Options n and m, it is preferable to perform as 
much of the processing as possible at the MR in order to reduce the amoimt of 
information transmitted from the MR to the Host. For example, if TDOA is 
employed for location, it is pref^able to measure the TOA of a received signal, 
ri(t), that is transmitted frorii the i*^ BST, at the MR The same is true for PDOA 
and FDOA. Ih the TQA case and in thp POA case, the requked processing 
consists of correlating the received signal with a clean version of Pi(t) and 
sending tihie information regarding the coirelatioii peaks to tiie Host In the JQA 
case, the required processing consists of correlatiag the received signal with 
several cleaii versions of the Pi(t) at different frequencies and sending the 
inforrnation regarding the correlation peaks to the Host 

The correlation peaks between ri(t) and a clean version of pi(t) can be super- 
resolved using SR algorithms. It is sometimes preferable to perfonn the SR 
algorithms at the Host rather than performing them at the MR since it requires 
ia number of computationally intensive operations that are easier to develop at 
flie Host than at the MR. This implies that the correlation peak needs to be 
transferred to the Host which does not represent a large amount of information 
and can be easily transferred to the Host. 

In order to be able to correlate the received signal, ri(t), at the MR, a clean 
version of Pi(t) must be available at the MR beforehand. In order to avoid 
correlating ri(t) with Pj(t), one can use either one of the following two 
procedures: 
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Procedure I: 

1 . Filter the entire received signal; 

2. Extract the phase of an existing tone. In this case, TOA and 
AOA information are contained in the phase of the tone while 
the FO A information are contained in the frequency of the tone 
and can be extracted using Procedure m (see below). 

Procedure 11: 

1 . Filter the entire received signal; 

2. Perform a nonlinear operation on the filtered signal; 

3. Extractthephaseof a tone generated by the nonlinear operation 
in step 2. In this case, TOA and AOA information are contained 
in the phase of the tone while the FOA information are 
contained in the frequenqr of the tone and can be extracted 
using Procedure DDL. 

The nonlinear operation in Procedure H can be either a quadratiG operation 
(e,g. squaring the. sigpal), a delay product operation (i.e. multiplying the si^^al 
by . a delayed version of itselQ or any other suitable noidmear operatioiL 

Procedures I and n cain be used both for acquiring the signal and for 
estimating its TOA. This is accomplished by stan^^ing the time using a reliable 
clock at either the start of the generated tone or obtaining the time stamp at the 
end of the observation window (the latter is preferred). 

Procedures I and EL can be improved by resolving the existing or generated 
tones using SR and inverse SR algorithms. Moreover, the filtering required in 
step 1 of Procedure n can be carried out in the digital domain in order to 
reduce the efifect of temperature and aging on group delay variation. 

Procedures I and II can be improved by repeating them a number of times and 
then either averaging the Location information or averaging the TDOA, PDOA 
or FDOA between pairs of BSTs. 

An alternative method for improving Procedures I and DC is by substantially 
increasing the duration of observation, Tj, of the received signal, ri(t), thereby 
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increasing the received SNR. In order to be able to increase T; at will, pj(t) has 
to be available at the MR. An important feature in Procedure II is the fact that 
the signal ri(t) transmitted by the i* BST does not have to be a deterministic 
signal, i.e. a clean replica of the signal does not have to be available at the MR. 
Equivalently, an important feature in Procedure II is the fact that the 
probability of location is 100%, unlike Procedure I which is below 100%. 
Hiis leads to the fact that in Procedure II, increasing the duration of 
observation, T,-, can acconunodate different situations. For example, in rural 
areas, where the cells are large, the SNR can be increased by increasing the 
duration of observation. Or equivalently, if the geometry is not good (i.e. 
HDOF is large, e.g. highway), one can improve it by increasing the duration of 
observation hence, increasing the number of BSTs that receive the signal ri(t) at 
a reasonable level. The duration of observation is thus preferably optimized 
with respect to the SNR and HDOP. In addition to reducing HDOP by 
increasing the number of BSTs, this also leads to reducing the effect of 
multipath, particularly in urba^ centers. 

In order to reduce the effect of in-band interference on the phase of the desired 
tone, narrowband digital filters of a few hertz bandwidth are used (e.g; Fourier 
Xransforitirbased filters, sliding Fourier Transform-based filters, IIR filters, 
FIR filters, etc.). The Fotirier Transfonn-based filter has a bandwidth which 
is a function of the observation duration, Tj. The IIR filter has a steady state 
response lime which is a function of Tj. In other words, the observation 
duration, Tj, has tp be selected in order to optimize the SNR (increase it) versus 
the bandwidth of the filter (decrease it) versus the response time of the wireless 
location system (minimize it). 

In order to reduce the effect of the multipath, MPi(t), on the performance of a 
kinematic wiretess location system, it is desirable to resolve the different 
Doppler frequencies Sfj (which correspond to different TOAs and dififerent 
AOAs) and to choose the frequency corresponding to the first arrival. The 
method of resolving the different frequencies can be accomplished using SR 
algorithms. 

Hardware 

In order to improve over prior art for modification II, one should attempt to 
1 . reduce or avoid: 
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— Gain Imbalance; 

— Phase Imbalance; 

— Carrier Feed-through; and 

— DC offset, 

2, be able to use the same RF front end for all standards (analog or 
digital) over a fixed band; and 

3. reject out-of-band Interferers (e.g. paging, trunked radio, etc.) using 
digital filters in the Digital Signal Processor (DSP). 

To achieve the three objectives above, it is proposed to use a linear IF- 
sampling receiver for the MR, or for the RR, or for both. Such an IF-sampling 
receiver is designed to have high rejection, low group delay variation and good 
sensitivity. In order to have high rejection with good sensitivity, the receiver is 
designed to have (see Figure 8): 
Design I: 

1. A number of Inteimediate Frequency (IF) stages (810) that convert the 
received RF sigpal Re{^(t)exp(j27tf;t)} (801) to an IF signal (807) of 
IF fceqaency f,. Initial rejection is achieved at RP using an RF B^id 
Pass Filter (BPF) (802), followed by an RF amplifier (803) for.good 
sensitivity. Subsequjent rejections are achieved at each IF stage which 
consist, of a itiixCT (804) followed by an Intamediate Frequency 0F) 
BPF filter (805) sind an IF amplifier (806). 

2. A pre-aiwsing filter (811) precedes the Analog-to-Digitai eonverter 
(812) that converts the IF analog signal to an IF digital signal. The 
sampling rate is intended to be below flie Nyquist rate in order to force 
aliased copies of the IF signal (807) close to a lower digital n^, fi- 

3. Digital filters (813) which reject noise and interference outside of the 
digital IF bandwidth. 

4. A Digital Down Converter (DDC) (804) which converts the (real) 
digital IF signal to a (complex) digital baseband signal (optional). 

In order to have low group delay variation, the following procedure is 

suggested: 

Procedure IH: 

1. Estimate the carrier offset (818), (due to Doppler, Sfj, and due to LO 
offset, Afi) of the signal ri(t) (801) using Digital Signal Processor 
(DSP) (817). 
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2. Adjust the Local oscillators (808) in the receiver corresponding to the 
MR during reception of rjCt) (801) based on the estimated carrier offset 
(818), in step 2 above. 

Ex planation of Procedure III: Since the overall group delay, Atj, is a function 
of frequency, it can vary significantly from one frequency to another. 
Furthermore, according to the AMPs standard, Afj can be as large as 2kHz. For 
this reason, it is important to force Atj to be independent of the carrier offset, 
A^- + 5^. This is accomplished by estimating the carrier offset, Afj + 5fj, in step 
1 of Procedure m and removing it in step 2 of Procedure m prior to 
processing it for location purposes. In other words, steps 1 and 2 force the 
overall group delay, Atj, to remain constant, thereby reducing overall group 
delay variations due to firequency offeets. In order to reduce overall group 
delay variations due to temperature and aging in the RR, a calibration 
procedure has to be performed on a regular basis on the RR- 

Notes; 

. 1. In Procedure m, step 1 can be performed using a discrete Carrier 
Recovery Loop (CRL) described in Procedure IV (see below). 

2, In i^rocedure jn, step 2 caa be performed using a Dir«:t Digital 
:^tli65iser 09DS) (809). . 

The discrete Carrier Recovery Loop consists of several steps that are repeated 
a number of times until a satisfectory level of pCTformance is reached: 
Procedure IV: 

1. Digitally Signal Process the sampled IF signal (814 or 816) usirig the 
DSP (817) in order to derive an adequate objective function, 4bj- 

2. If is optimized to within a certain constraint, stop the loop, 
otherwise: 

3. search for a new freqiienQr offset (818), 

4. adjust the LO (808) by the new frequenqr offset (818), and 

5. go back to step 1. 

In Procedure IV, step 3 can be performed using a gradient-type search 
algorithm or any other suitable search algorithm. 
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In Procedure IV, step 4 can be performed using either a DDS (809) only or a 
DDS and a Digital Down-Converter (DDC) (815). If a DDS (815) only is used 
in step 4, then step 1 in Procedure IV processes a new segment of the signal, 
ri(t), every iteration. If a DDS (809) and a DDC (809) are used is step 4, then 
step 1 in Procedure IV processes the same stored segment of the signal, ri(t), 
every iteration until the DDS (809) is used. It is preferred to use both the DDC 
and the DDS in order to reduce the effect of finite precision arithmetic on the 
estimation process. 

In the case when Design I does not contain step 4, the DSP (817) in step 1 of 
Procedure IV of the carrier recovery loop can consist of a digital notch filter 
centered at the desired IF frequency f^, whose power of response is to be 
minimized, i.e. i^j is the response of the digital notch filter. In this case, the 
notch/filter can consist of 

1. a first digital Band Pass Filter centered at the fj H- firequency of 
desired Tone; and 

2. a second digital Band Pass Filter centered at the - fi?equency 
of desir^ Tone. 

The Band Pass filtqr$ can be iixQ>lemented using , (real} digital Biquadratic 
filters. 

Alternatively^ the DSP (817) in step 1 in Procedure IV of the carrier 
recovery loop can consist of a digital Band Pass Filter centered at the desired 
IF fi:cquency fz, whose power oif response is to be maximized, i.e. is the 
response of the Band Pass Filter. In this case, the BPF can be implemented as a 
(coniplex) 1^ order (or higher) digital IIR filter. 

In the case when Design I contains step 4, the signal to process is a baseband 
signal (816) and the digital BPF tums into a simple (complex) LPF. An 
example of such a LPF is an averaging filter which adds all real samples 
together and all imaginary samples together over the observation interval. 

In order to reduce the effect of clock errors and LO of&et (or drift), A^, a 
conmion reference is used for all the LOs and for all the clocks. A suitable 
common reference for the LOs and for the clocks can be derived fi*om the GPS 
signal. 
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Given that either all LOs in the BSTs have a common reference, which implies 
that 

Af, = Af2 = Af3 = Af 

in equations (32), or that the RR estimates Af„ Afj, Af3 in order for them to be 
removed, it is possible to estimate the speed, v, with respect to the BSTs and 
Direction Of Travel (DOT), <p, of the MR from the carrier offsets that are 
estimated from Procedure V at each BS using the discrete carrier recovery 
loop. This can be accomplished using the Frequency Difference of Arrival 
(PDOA) between each pair of BS as expressed in eqiiations (33). 

Calibration methods at the RR 

It is preferred to ininiinize the use of the RRs in order to reduce cost and 
complexity of a handset-based WLS, It is desired to produce a handset-based 
Wireless Location System (WLS) where an MR can be located without the 
need for an RR- This is achieved by considering a transmitter in an existing 
cellular sector as . a BST. Hie advantage is that such a sector is usually time- 
refcrraced to other sectors belon^g to flie same .SST, i.e. sectors in the same 
. BS are iriore or less synchronized to one another. This implies that the MR can 
estimate the AOA of sigcials transmitted by such sectors without a need for an 
RR. Moreover, the MR can estimate its location using only two independent 
BSs based on the individual estiniated AOA froni each BS withdut the need for 
anRR- 

A wireless location system has to be calibrated prior to operation (and 
occasionally during operation) in order to remove the effect of the overall 
Group Delay, Atj, at each RR- Moreover, temperature and aging can cause the 
overall Group Delay to change with time: Therefore, it is necessary to calibrate 
the RRon a regular basis. A method to calibrate an RR is as follows: 
Procedure V: 

1. transmitting at precise time intervals signals from a 
transmitter of known location to the RR over several frequency 
bands. Both the RR and the transmitter have to be referenced to 
the same clock source; 

2. characterizing the frequency response of the RR; 
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3. estimating the time of travel between the transmitter and the 



RR; 

4. estimating the group delay of the RR. 



Procedure V can be repeated a number of times in order to average the relative 
group delay over time. In order to reduce deployment cost of a transmitter at 
known locations, an MR is stored in proximity to the RR. This allows for 
automatic calibration without a need for human intervention by forcing the MR 
. to act as a transmitter. 

In order to minimize the relative group delay variations between calibration 
intervals, it is recommended to use SAW filters as die BPF (803, 805, 81 1) in 
Figure 8. It is also recommended to place the SAW filters (803, 805, 81 1) in a 
metallic block to mininiize flexing effects due to temperature variations. 

HostjSoftware 

CHven a nxmiber of positioning algorithms and a nxunber of methods for 
blunder detection, positioning strategy is used as Ulustrated in Figures 9a, 9b, 
9c and 9d. Figures 9 illustrate ffie preferred embodiriient fof adjusting the 
observed TOAs- such that the sum of the squares of the residuals is a minimum. 
Other mefliods for adjusting the observed according to some optimization 
CTiteria aire available including: 

• General Minimimi Variance Unbiased Estimation, 

• Best Linear Unbiased Estimation 

• Maximimi Likelihood Estimation, 

• Method of Moments, 

• General Bayesian Estimation, 

• Linear Bayesian Estimation, 

• Kalman Filtering, etc. 

The positioning process begins with a set of TOA observations from a 
minimum of tiiree BSs. If only three TOAs are available (901), solution 
bifiu-cation is tested for (902). In the event of solution biftircation, the LOCA 
algorithm (903) is performed on the TOA observations yielding two solutions. 
These two solutions are each used as the initial position for Least Squares (904 
and 905) resulting in two Least Squares solutions for the MR. The Least 
Squares algorithm here, and in all other parts of Figure 9, is Hybrid TDOA 
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positioning with received signal strength used to weight the individual BS 
observations. During every iteration of Least Squares, misclosures for each 
BST are calculated. If any misclosure exceeds some multiple of the RMS 
value of the entire set of misclosures, the corresponding BST, and its TOA, is 
no longer used within that particular execution of Least Squares. 

If both Least Squares solutions (904 and 905) converge (906) but are not equal 
(907), it is assumed that Least Squares has converged to the two possible 
solutions and both Least Squares solutions are reported as possible positions 
for the MR (908). If the either of the LS solutions diverge (906) or the 
convergent solutions are equal (907), the two LOCA solutions from (903) are 
reported as possible positions of the MR (909). 

If at (902) bifurcation is found to not exist, LOCA is again executed (910) but 
yields only one solution. This solution is used as the initial position in Least 
Squares (911). 

When redundancy does exist (901), Plane hitersection (912) (or any other 
closed-form position estimation algorithm) is executed. The resulting solution 
is used as the initial position for Least Squares (913). At this point a check is 
made as to whether Least SquareS: (913 or 911) converged (914). If Least 
Squares converges, another check for observational redundancy is made (915), 
If there is no redimdancyi the Least Squares solution is reported as the position 
estimate of flie MR (1923). If fliere is redundancy, the normalized rbsiduaU arfe 
statistically tested for normality (916). The Least Squares solution is reported 
(923) should all residuals pass. If any residuals fail, the redimdancy numbers 
of the failing observations are checked (917). The redundancy number of the i*^ 
observation is defined as 

g,=(cA-')i, (52) 
where Cy is the covariance matrix of the residuals and is defined as 

C,=C,-a[a^C,-'a]''a^. (53) 
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Should all failing observations have redundancy numbers less than some 
threshold (preferably, this threshold is 0.5), the Least Squares solution is 
reported as the MR position estimate (923). 

If observations with failing residuals have redundancy numbers larger than the 
threshold (917), that observation with redundancy number greater than the 
threshold and with the largest standardized residual is permanently removed 
from the observation set (928). The initial position used for the previous 
execution of Least Squares is remembered (929) and used again in Least 
Squares with the truncated observation data set (930 or 93 1). 

Should Least Squares now diverge (932), the previous Least Squares solution 
which did converge is reported as the position estimate of the MR (936). If 
Least Squares does converge (932), and there is no redundancy (933), the 
newly convergent Least Squares solution is reported (936). If there is 
redundancy (933), the standardized residuals are tested for normality (934). 
Should all the residuals pass, the newly convergent Least Squares solution is 
reported (936). Otherwise, if all failing observations have redundancy numbers 
less than some threshold preferably, this tiireshold is 0.5), flie newly 
convergent Least Squares solutioti is reported (936). 

If observations with failing residuals have redundancy numbers largor tiian the 
threshold (935), that observation with redundancy nmnber greater than the 
threshold and witii the largest standardized residuail is permanenfly removed 
from the observation set (928). The process then continues as described 
iminediately above. 

If at (914) Least Squares diverges. Least Squares is executed again but with the 
average of participating BSTs coordinates used as the initial position (918). If 
Least Squares now converges (919) but there is no observation redundancy 
(920), the newly convergent Least Squares solution is r^orted as the position 
estimate of the MR (923). If there is redundancy (920), tiie standardized 
residuals are tested for normality (921). Should all residuals pass, the newly 
convergent Least Squares solution is reported (923). If some of the residuals 
fail but none of the corresponding observations have redundancy numbers 
greater than some threshold (922), the newly convergent Least Squares 
solution is reported (923). 
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If observations with failing residuals have redundancy numbers larger than the 
threshold (922), that observation with redundancy number greater tiian ttie 
threshold and with the largest standardized residual is permanently removed 
from the observation set (928). The process then continues from (928) as 
described above. 

When Least Squares does not converge at (919), and there is no redundancy 

(924) , either no solution is reported for this particular set of observation data or 
the solution from the last iteration before divergence is reported (927). If there 
is redundancy, but all of the standardized residuals pass the normality test 

(925) , no solution is reported for this particular set of observation data or the 
solution from the last iteration before divergence is reported (927). Should 
some of the residuals fail, the obiservation with the largest standardized residual 
is permanently removed from the data set regardless of its redundancy number. 
The process thra begins at point (912) as described above. 

Figures 9 . are also applicable for, 

• AOApositiomng, 

• AOA/TnDOApojsitioning, and 

• AOA/Range positioning. 

Transferring I/Ocatioh Information of the MR over th e Internet 
Given that the Intemet is global and inexpensive, the communication between 
the Host and the customer can be achieved over it For example when a MR 
calls 911 for emergency, it is possible to relay the positional information of the 
MR from the Host to the PS AP over the intemet. Similarly, in the case when a 
powered-on MR is to be located by a customer, its positional information can 
be relayed from the Host to the customer through the intemet. 

The utility of the intemet in its capacity of transferring MR location 
information on a global scale, extends the application of the wireless location 
system beyond the cellular network for which the BSTs have been deployed. 
Through the latest intemet technology such as Java, JavaBeans, as well as 
CORBA (Common Object Request Broker Architecture), MR location 
information residing at the Host can be integrated with third-party infomiation 
(i.e. a map database, or a database consisting of geo-coded business addresses 
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such as restaurants, towing companies, etc.) residing in some other geographical 
location, perhaps even in another country. The combination of this MR location 
information with the third party information can be treuisferred over the internet 
to allow customers to locate a particular MR wrt to either a map location (i.e. 
street address) and/or a business location, provided that the MR is within 
network coverage for the wireless location system. Through this process the 
MR location information and the third party database information can be 
accessed or "served up" to the customer through eflScient Java Internet 
technology processes. The service provider will integrate together the various 
components including the MR location information derived from the wireless 
location system and the third-party information database. The integration process 
of the different databases is transparent to the customer. The customer will only 
know that the location-based service exists as such to bring all of the different 
components together to provide for a complete service which can be offered 
either on a regional, national or perhaps global basis. Through the internet these 
location-based services can now be offered on an economical basis to the 
customer. Examples of such services are fleet management, concierge services, 
roadside assistance, child find services, etc. Legality and security aspects are a 
concern on the internet and in some cases a dedicated link may be necessary. 

A clean replica refened to herein is a copy of a signal that is not noisy, 
as for example may bq obtained by a stored version of the sigoal transmitted. 
Low group delay variation preferably means a few 10s of nanoseconds of group 
delay. An adequate objective fiinction is a fimction known in the literature, for 
example minimizing the mean squared error. A PSAP is a public safety 
answering point, for example a 91 1 station. An RSSI is a receiver signal istrength 
indicator. 

Every document referred to herein as containing useful information for 
implementing MR location is incorporated herein by reference. 

Modifications to the preferred embodiments disclosed may occur to 
those skilled in the art without departing from the scope of the invention. 
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We Claim: 

I . A method of estimating the location of a Mobile Receiver (MR) comprising 

the steps of: 

5 receiving a plurality of signals from a plurality of Base Station Transmitters 

(BSTs) by the MR; and 

performing a nonlinear operation on the plurality of received signals at the 
MR in order to generate at least one tone per received signal. 

10 2. The method of claim 1 wherein the nonlinear operation is selected from the 

group consisting of a quadratic operation, a delay product operation; and an absolute value 
operation. 

3, The method of claim 1 where the location operation is performed by 
1 5 estimating the phases of the generated at least one Tone. 

4. The method of claim 1 wherein the received signals are filtered prior to 
performing the nonlinear operation. 

20 5. The method of claim 4 wherein the received signals are correlated with a 

clean replica of themselves prior to,performing the filtering operation. 

6. A method of estimating the location of a Mobile Receiver (MR) comprising 

the steps of: 

25 receiving a plurality of signals from a plurality of Base Station Transmitters 

(BSTs) by the MR; and 

estimating the location of the MR by estimating the phases of existing tones 
in each received signal. 

30 7. The method of claims 1 or 6 where the generated or existing Tones are 

resolved using SR or Inverse SR algorithms. 

8. The method of estimating the location of the MR in claims 1 or 6 by 

extracting the Tones using digital filtering. 

35 
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9. The method of claim 8 wherein the method of digital filtering comprises the 

steps of 

performing a Fourier Transform (FT) on segments of the signal, 
windowing the desired frequency domain bands for each FF segment; and 
performing an inverse Fourier Transform (IFT) on the windowed bands. 



10. The method of claim 9 wherein the windowed frequency domain bands for 

each FT segment are transferred to a Host. 

10 11. The method of claim 8 wherein the method of digital filtering performs an 

operation selected from the group consisting of:a sliding (or recursive) FT operation, 
an IIR filtering operation; and 
an FIR filtering operation. 



15 12. The method of claim 10 where a time domain signal, ri(t), is constructed at 

the Host by 

appending all the received signal frequency domain segments; and 
performing an IFT operation. 

20 13. The method of claims 9 or 12 where a plurality of tones in the reconstructed 

time domain signal, ri(t), are resolved using SR or Inverse SR algorithms. 



14. The method of claim 13 where a tone is selected from the plurality of 

resolved tones such that it corresponds to the earliest Time-Of-Arrival among all resolved 
25 tones. 



15. The methods of claims 7 or 13 where the SR or the inverse SR algorithms is 
selected from the group of algorithms consisting of: MUSIC, ESPRIT, Auto Regressive 
Moving Average, Minimum Variance, MUSIC using Higher Order Statistics, ESPRIT 

30 using Higher Order Statistics, Auto Regressive Moving Average using Higher Order 
Statistics, and Minimum Variance using Higher Order Statistics. 

16. The method of claim 1 5 further comprising 

the method performing an Inverse SR algorithm; 
35 the inverse SR algorithm comprising: 

a conventional lime domain correlator. 
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a lime domain window, 
a Fourier Transform, 
a frequency domain window, 
a frequency domain equalizer, and 
5 a processor which performs a SR algorithm in order to resolve the TOAs 

in the received signal at a given BST. 

1 7. The method of claim 1 5, further comprising 

the method performing an Inverse SR algorithm; 
10 the inverse SR algorithm comprising: 

a frequency domain correlator, 
a frequency domain window, 
an inverse Fourier Transform, 
a time domain window, 
1 5 a time domain equalizer, and 

a processor which performs a SR algorithm in order to resolve the FOAs 
in the received signal at a given BST. 

20 1 8. A method of estimating the location of a MR from a plurality of signals 

transmitted from a plurality of BSTs, where one or more antennas is used at the MR and in 
which location of the MR is achieved using a technique selected from the group consisting 
of:: 

AOA positioning (multi-angulation), 

25 TOA positioning (circular multilateration). 

Range positioning (circular multilateration), 
TDOA positioning (hyperbolic multilateration), 
TDOA positioning (circular multilateration with TDOAs), 
AOA/TDOA positioning (multi-angulation/hyperbolic multilateration), 

30 AOA/TDOA positioning (multi-angulation/circular multilateration with 

TDOAs), 

AOA/Range positioning (multi-angulation/circular multilateration), 
AOA/Range positioning (multi-angulation/circular multilateration), and 
combining the signals from all antennas for location purposes. 

35 

19. The method of claim 18 in which more than one antenna is used. 
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20. The method of eslimating ihc location of a MR in claim 18 where the 
combining of the signals from all antennas is selected from the group consisting of: 

selection combining, 
5 maximal ratio combining, 

co-phasing combining, and 
equal gain combining. 

21. A method of estimating the location of a MR from a plurality of signals 
10 transmitted from a plurality of BSTs, using equation (29) and equation (30). 

22. The method of claim 18 where location of the MR is achieved using a 
technique selected from the group consisting of: 

Range positioning (circular multilateration) with Round Trip Delay, and 
15 AO A/Range positioning (multi-angulation/circular multilateration) with 

Round Trip Delay. 

23. A method of estimating the location of a MR wherein location of the MR 
is achieved based on a solution selected firom the group consisting of: 

20 solution I, 

solution II, 
solution III, and 
solution IV. 

25 24. A system to estimate the location of a Mobile Receiver (MR) comprising 

of: 

a receiver at the MR for receiving a plurality of signals from a plurality of 
Base Station Transmitters (BSTs); and 

a computing means connected to the receiver by a communication linkfor 
30 performing a nonlinear operation on the plurality of received signals in order to generate at 
least one tone per received signal. 

25. A system to estimate the location of a Mobile Receiver (MR) comprising of: 

a receiver at the MR for receiving a plurality of signals from a plurality of 
35 Base Station Transmitters (BSTs); and 
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a computing means connected by a communication link to the receiver for 
estimating the location the MR by estimating the phases of existing tones in each received 
signal. 

5 26. The system of claim 24 further comprising a computing means for 

performing an Inverse SR algorithm; 

the inverse SR algorithm comprising: 

a conventional time domain correlator, 
a time domain window, 
10 a Fourier Transform, 

a frequency domain window, 

a frequency domain equalizer, and 

a processor which performs a SR algorithm in order to resolve the TOAs 
in the received signal at a given BST. 

15 

27. The system of claim 24 further comprising a computing means for 
performing an Inverse SR algorithm; 

the inverse SR algorithm comprising: 

a firequency domain correlator, 
20 a firequency domain window, 

an inverse Fourier Transform, 
a time domain window, 
a time domain equalizer, and 

a processor which performs a SR algorithm in order to resolve the FOAs 
25 in the received signal at a given BST. 

28. A system to estimate the location of a MR from a plurality of signals 
transmitted from a plurality of BSTs, where one or more antennas are used at.the MR and 
wherein the system comprises means for implementing a technique selected from the group 

30 consisting of: 

AOA positioning (multi-angulation), 
TOA positioning (circular multilateration). 

Range positioning (circular multilateration), 
TDOA positioning (hyperbolic niullilateralion), 
35 TDOA positioning (circular multilateration with TDOAs), 

AOA/TDOA positioning (multi-angulalion/liyperbolic multilateration). 
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AOaAI^DOA positioning (miilli-angulation/circuiar mullilateralion with 

TDOAs), 

AOA/Range positioning (muUi-anguIalion/circular mullilateralion), 
AOA/Range positioning (multi-angulaiion/circular multilateration), and 
5 combining the signals from all antennas for location purposes. 

29. The system of claim 28 in which more than one antenna is used at the 

MR. 

10 30. A system to estimate the location of a MR from a plurality of signals 

transmitted from a plurality of BSTs, using equation (29) and equation (30). 

31. The system of claim 28 further comprising means for implementing a 
technique selected from the group consisting of: 

15 Range positioning (circular multilateration) with Round Trip Delay, and 

AOA/Range positioning (multi-angulation/circular multilateration) with 
Round Trip Delay. 

32. A system to estimate the location of a MR from a plurality of signals 
20 transmitted from a plurality of BSTs, based on any one of the following solutions: 

solution I, 
solution II, 
solution in, or 
solution IV. 

25 

33. A system to estimate the location of a Mobile Receiver (MR), the system 
comprising a Receiver that receives a plurality of signals transmitted from a plurality of 
BSTs; 

30 wherein the Receiver includes an IF-sampling Receiver comprising: 

a number of IF stages to convert a received RF signal n(t) to an 
analog I F signal of IF frequency f i ; 

a pre-aliasing filter; 

an Analog-to-Digital (A/D) converter thai converts the IF analog 
35 .sigiicil lo an IF digital signal; and 

a processor for processing the IF digital signal. 
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34. The system of claim 33 in which each IF stage of the IF-sampIing receiver 
comprises: 

a mixer; 

5 an IF Bandpass Filter (BPF) filter; and 

an IF amplifier. 

35. The system of claim 33 wherein the sampling rate of the AJD converter is 
below the Nyquist rate in order to force aliased copies of the digital IF signal close to a 

1 0 lower Intermediate Frequency f2- 

36. The system of claim 33 wherein the processor performs functions selected 
from the group consisting of the functions: 

digital filtering to reject noise and interference outside of the digital IF 

15 bandwidth; 

digital down conversion to convert the digital IF signal to a complex digital 
baseband signal; and 

discrete carrier recovery in a carrier recovery loop. 

20 37. The system of claim 33 wherein the IF-sampling Receiver has low group 

delay variations over a desired frequency band. 

38. The system of claim 33 where the IF-sampling Receiver comprises one or 
more local oscillators (LOs) and clocks and has a common reference for all the LOs and for 

25 all the clocks. 

39. The system of claim 38 wherein the common reference for the LOs and for 
the clocks is GPS-based, 

30 40. The system of claim 34 wherein the pre-aliasing filter is a SAW filter and 

the BPF is a SAW filter. 

41. The system of claim 40 wherein the SAW filters are placed in a metallic 
block to minimize flexing or temperature effects. 

35 

42. The system of claim 33 in which the Receiver is the MR. 
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43. The system of claim 33 in which the Receiver is a Reference Receiver 

distinct from the MR. 

5 44. A system to estimate the location of a MR using a Reference Receiver 

having a LO, the system comprising: 

an estimator for estimating the carrier offset (due to Doppler, Sfi, and due to 
LO offset, Afj) of the plurality of signals transmitted by a plurality of BSTs and received by 
the RR, 

10 

45. ■ The system of claim 44 flirther comprising a discrete Carrier Recoveiy 
Loop. 

46. The system of claim 45 wherein the carrier recovery loop comprises a 
15 processor configured to repeat the following steps until a satisfactory level of performance 

is reached: 

1. process the sampled IF to derive an adequate objective function, ^bj; 

2. if £,b is optimized to within a certain constraint^ stop the loop, otherwise: 

3. search for a new frequency offset, 
20 4. adjust the LO by the new frequency offset, and 

5. go back to step 1. 

47. The system of claim 46 wherein the search for a new frequency offset in 
step 3 is performed using a gradient-type search algorithm. 

25 

48. The system of claim 46 wherein adjusting the LO by the new frequency 
offset can be performed either using a DDS or a Digital Down-Converter or a combination 
thereof. 

30 49. The system of claim 36 or 46 wherein the processor performs the function 

of a carrier recovery loop by perfonning a digital notch filtering operation centered at the 
desired IF frequency fi, and whose power of response is an optimized objective function, 
fohjr derived from the sampled IF. 

35 50. The system of claim 49 in which the digital notch filtering operation 

consists of: 
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a first digital Band Pass Filter centered ai the F2 + frequency of desired 

Tone; and 

a second digital Band Pass Filter centered at the f2 - frequency of desired 

Tone. 

5 

51. The system of claim 36 or 46 wherein the processor is configured to 
perform the function of a carrier recovery loop by performing a digital Band Pass filtering 
operation centered at the desired IF frequency f2, and whose power of response is an 
optimized objective function, £>bj, derived from the sampled IF signal. 

10 

52. . The system of claim 51 wherein the desired IF frequency f2 is equal to zero 
and the digital BPF is a LPF. 

53. The system of claim 51 wherein both Band Pass filters are implemented 
1 5 using digital B iquadratic filters. 

54. The system of claim 46 further comprising a Host and wherein the estimated 
carrier offset last found is transferred to the Host in order to aid in estimating the speed v of 
the MR and its Direction Of Travel. 

2Q 

55. A system to estimate the location of a MR having a LO, the system 
comprising: 

an estimator for estimating the carrier offset (due to Doppler, Sfi, and due to 
LO offset, Afi) of a plurality of signals transmitted by a plurality of BSTs and received by 
25 the MR. 

56. The system of claim 55 in which the estimator comprises a discrete Carrier 
Recovery Loop. 

30 57. The system of claim 56 wherein the carrier recovery loop comprises a 

processor configured to perform the following steps until a satisfactory level of 
performance is reached: 

1 . process the sampled IF to derive an adequate objective function, fobj; 

2. if fob is optimized to within a certain constraint, stop the loop, otherwise: 
35 3. search for a new frequency offset, 

4. adjust the LO by the new frequency offset, and 
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5. go back lo step I. 

58. The system of claim 57 wherein the search for a new frequency offset in 
step 3 is performed using a gradient-type search algorithm. 

5 

59, The system of claim 57 wherein adjusting the LO by the new frequency 
offset is performed either using a DDS or a Digital Down-Converter or a combination 
thereof. 

10 60. The system of claim 46 wherein the processor is configured to perform the 

function of a. carrier recovery loop by performing a digital notch filtering operation 
centered at the desired IF frequency fi, and whose power of response is an optimized 
objective function, fobj, derived from the sampled If signal. 

15 61. The system of claim 59 in which the digital notch filtering operation 

consists of: 

a first digital Band Pass Filter centered at the fz + frequency of desired 

Tone; and 

a second digital Band Pass Filter centered at the fa - frequency of desired 

20 Tone. 

62. The system of claim 61 wherein the processor is configured to perform the 
function of a carrier recovery loop by performing a digital Band Pass filtering operation 
centered at the desired IF frequency fa, and whose power of response is an optimized 

25 objective function, fobj, derived from the sampled If signal. 

63. The system of claim 62 wherein the desired IF frequency (2 is equal to zero 
and the digital BPF is a LPF. 

30 64, The system of claim 61 wherein both Band Pass filters are implemented 

using digital Biquadratic filters. 

65. The system of claim 55 or 56 further comprising a Host and wherein the 

estimated carrier offset last found is transferred to the Host in order lo estimate the speed v 
35 of the MR and its Direction Of Travel. 
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66. The system of claim 65 wherein ihe speed of ihe CT and its direction of 
traveling is estimated as a function of the difference between the estimated carrier offsets. 

67. A method of estimating the location of a MR, the method comprising the 
5 steps of: 

receiving a plurality of signals transmitted from a plurality of BSTs at a Receiver; 
converting the received RF signal ri(l) to an IF signal of IF frequency f| in an IF- 
sampling Receiver; 

filtering the IF signal with a pre-aliasing filter; 
10 converting the filtered IF analog signal to an IF digital signal in an Analog- 

to-Digital converter;, and 

processing the IF digital signal in a signal processor. 

68. The method of claim 67 in which the Receiver is the MR. 

15 

69. A method of estimating the location of a MR, the method comprising the 
steps of; 

estimating the carrier offset (due to Doppler, 5fi, and due to LO offset, Afi) 
of a plurality of signals transmitted by a plurality of BSTs and received by a RR. 
20 . 

70. A method of estimating the location of a MR, the method comprising the 
steps of: 

estimating the carrier offset (due to Doppier, Sfi, and due to LO offset, Afj) 
of a plurality of signals transmitted by a plurality of BSTs and received by the MR. 

25 

71. A method of estimating the location of a MR from a plurality of signals 
transmitted from a plurality of BSTs, the method comprising the steps of: 

characterizing the relative group delay, Ati - Atj, at a RR between a signal 
received from the i**" EST and a signal received from the j'^ BST prior to operation; 
30 removing the characterized relative group delay, Ati - Atj, from the 

estimated Time Difference Of Arrival, x-, - xj, at the MR between a signal received from the 
i*** BST and a signal received from the j'** BST to yield a corrected Time Difference of 
Anival; and 

Imdin«i the location of the MR using the corrected Time Difference of 

35 Anival. 
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72. The method of estimating the location of a MR in claim 71 wherein 

characterizing the relative group delay, At; - Alj is achieved by 

forcing a transmitter of known location to transmit at a known time of 
transmission (TOT); 
S receiving the transmitted signal at the RR; 

measuring the time of arrival (TOA) of the signal at the RR; 

estimating the time of travel as the difference between the TOA and the 

TOT; and 

finding the difference between the estimated time of travel and the known 
10 time of travel, the difference being the estimated delay of the transmitted signal through the 
RR. . 



73. The method of claim 73 further comprising repeating the steps of claim 72 
a plurality of times to yield a plurality of estimated delays, then averaging the estimated 

15 delays. 

74. The method of claim 72 wherein the transmitter of known location is a 
CT. 

20 75. A method of estimating the location of a MR in location state 0 or location 

state 1 from a plurality of signals transmitted from a plurality of BSTs, wherein the 
transfer of the positional information from the Host to a user is achieved over the internet. 

76. The method of claim 75 in which the MR is calling 91 1 (local emergency 
25 telephone number) and the information is transferred to a PSAP. 

77. The method of claim 75 further comprising extending the capability of a 
handset-based wireless location system beyond the network into which it is .deployed by 
using an internet technology as the basis to efficiently integrate third-party database 

30 information from various locations to provide holistic location-based services on a 
regional, national or global basis. 

78. The method of claim 77 wherein the internet technology is Java based. 

35 79. A method of estimating the location of a MR from a plurality of signals transmitted 
from a plurality of BSTs, comprising: 
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using Least Squares to solve for the positional information of the MR and 
using TDOA observations from distant BSTs to resolve for ambiguities in an AOA 
solution; 

using AOA observations from different sectors at the same cellular BS to resolve for 
S ambiguities in an AOA solution; and 

using range observations from different BSTs to resolve for ambiguities in an AOA 
solution. 

80. The method of claim 79 further comprising quantifying the effect of 
10 geometry on position estimation, with both AOA and range observations, using the design 

matrix in equation (48). 

81. A system for estimating the location of a MR from a plurality of signals 
transmitted from a plurality of BSTs, the system including a host, and the system 

15 comprising; 

the MR being adapted to receive some or all the following information from a 
remote source and to apply said information to the execution of a location state: 

identification of base station pilot signals for the mobile to search for pilot 
detection threshold; 

20 duration of time for searcher to accumulate pilot energy at each delay offset; 

search window time length or duration; 

spacing of successive sampling points within the search window; 
frequency at which to conduct the search for pilot signals; 
cellular system which the mobile should acquire; 
25 the number of times the mobile should repeat all or part of the location 

procedure; 

correlation interval; 
locations of base stations; and 

the MR being adapted to transmit one or more of the following information 
30 to the Host: 

times of arrival of pilot signals; 

identification of the pilot signals reported by said device; 
received signal qualities of pilot signals; 
pilot search results; and 
35 carrier frequencies of reported pilot signals. 
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82. A system for estimating the location of a MR from a plurality of signals 

transmitted from a plurality of BSTs, the system comprising: 

first computing means for characterizing the relative group delay, At; - Atj, at 
a RR between a signal received from the i'** BST and a signal received from the j**' BST 
5 prior to operation; and 

second computing means for removing the characterized relative group 
delay» Ati - Atj, from the the estimated Time Difference Of Arrival, x, - x,, at the MR 
between a signal received from the i*** BST and a signal received from the BST. 

10 83. A system for estimating the location of a MR in location state 0 or in 

location slate 1 from a plurality of signals transmitted from a plurality of BSTs, wherein 
the transfer of the positional information from the Host to a user is achieved over the 
internet. 

15 84. The system of claim 83 in which positioned data from a Host of a MR 

calling 91 1 is transferred over the internet. 

85. A system for estimating the location of a MR from a plurality of signals transmitted 
from a plurality of BSTs, where: 

20 Least Squares is used to solve for the positional information of the MR; 

TDOA observations from distant BSTs are used to resolve for ambiguities in an 
AOA solution; 

AOA observations from different sectors at the same cellular BS are used to resolve 
for ambiguities in an AOA solution; and 
25 AOA observations from different sectors at the same cellular BS are used to resolve 

for ambiguities in an AOA solution. 

86. The system of claim 85 further comprising a first computing means for 
quantifying the effect of geometry on position estimation, with both AOA and range 

30 observations, using the design matrix in equation (48). 

87. A method to estimate the location of a MR from a plurality of signals 
transmitted from a plurality of BSTs; 

the MR receiving some or all the following information from a remote source and 
35 applies said information to the execution of a location state: 
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idemificalion of base station pilot signals for the mobile lo search for pilot 
detection threshold; 

duration of time for searcher to accumulate pilot energy at each delay offset; 
search window time length or duration; 
5 spacing of successive sampling points within the search window; 

frequency at which to conduct the search for pilot signals; 
cellular system which the mobile should acquire; 

the number of times the mobile should repeat all or part of the location 

procedure; 
10 correlation interval; 

locations of base stations; and 
wherein some or all of the following information is transmitted from the MR to the 

Host: 

times of arrival of pilot signals; 
15 identification of the pilot signals reported by said device; 

received signal qualities of pilot signals; 
pilot search results; 

carrier frequencies of reported pilot signals. 

20 88. A method of estimating the location of a MR from a plurality of signals 

transmitted from a plurality of BSTs, where misclosures and/or standardized residuals are 
used to flag observations that might contain a blunder. 

89. The method of claim 88 in which location information from each BST is 

25 weighted based on its corresponding RSSI and wherein a method selected from the group 
consisting of the following methods is used to solve for the positional information of the 
MR: 

General Minimum Variance Unbiased Estimation, 
Best Linear Unbiased Estimation 
30 Maximum Likelihood Estimation, 

Least Squares Estimation, 

Method of Moments, 
General Bayesian Estimation, 
Linear Bayesian Estimation, and 
35 Kalman Filtering. 
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90. The method of claim 89 wherein a method sclecled from the group 
consisting of the following methods is used to solve for the positional information of the 

CT: 

Genera] Minimum Variance Unbiased Estimation, 
5 Best Linear Unbiased Estimation 

Maximum Likelihood Estimation, 
Least Squares Estimation, 
Method of Moments, 
General Bayesian Estimation, 
i 0 Linear Bayesian Estimation, and 

Kalman Filtering* 

91. The method of claim 89 wherein the flagged observations are either 
discarded or kept in the LS solution based on their effect on the total residuals of the LS 

15 solution and their individual redundancy number. 

92. The method of claim 91 where Least Squares is repeated without the 
flagged observations. 

20 93. The method of claim 92 where newly flagged observations are removed 

and Least Squares repeated until Least Squares converges or no more redundant 
observations exist. 

94. The method of claim 88 which uses a method selected from the group 
25 consiting of Chaffee's method and LOCA to detect solution bifurcation. 

95. The method of claim 93 where two estimates of the MR position are 
provided in the event that solution bifurcation does exist and there is no observational 
redundancy. 

30 

96. The method of claim 88 in which the hybrid TDOA positioning model in 
equation (42), which yields misclosures and residuals for the individual BSTs, is used. 

97. The method of claim 96 wherein a method selected from the group 
35 consisting of the following methods is used to solve for the positional information of the 

MR; 
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General Minimum Variance Unbiased Estimalion, 
Best Linear Unbiased Estimation 
Maximum Likelihood Estimalion. 
Least Squares Estimation, 
5 Method of Moments, 

General Bayesian Estimation, 
Linear Bayesian Estimation, and 
Kalman Filtering. 

10 98. A method of estimating the location of a MR from a plurality of signals 

transmitted from a plurality of BSTs, where a combination of AOA and TDOA 
observations are used such that the MR can be located with as few as two AOA 
independent observations and one TDOA observation. 

99. The method of claim 98 wherein a method selected from the group 
consisting of the following methods is used to solve for the positional information of the 
MR: 

General Minimum Variance Unbiased Estimation, 
Best Linear Unbiased Estimation 
Maximum Likelihood Estimation, 
Least Squares Estimation, 
Method of Moments, 
General Bayesian Estimation, 
Linear Bayesian Estimation, and 
Kalman Filtering. 

100. A method of estimating the location of a MR from a plurality of signals 
transmitted from a plurality of BSTs, wherein a combination of AOA and TOA 
observations are used such that the MR can be located with as little as one AOA 

30 observation and one TOA observation. 

101. The method of claim 100 wherein a method selected from the group 
consisting of the following methods is used to solve for the positional information of the 
MR: 

35 General Minimum Variance Unbiased Estimation, 

Best Linear Unbiased Estimation 



20 
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Maximum Likelihood Estimation, 
Least Squares Estimation, 
Method of Moments, 
General Bayesian Estimation, 
5 Linear Bayesian Estimation, and 

Kalman Filtering. 

102. The method of claim 100 where the TOA observations are based on 
estimating the RTD between the BS transmission and the reception of the MR*s response 

10 to the BS transmission. 

103. A method of estimating the location of a MR froni a plurality of signals 
transmitted from a plurality of BSTs, where a combination of TOA observations are used 
such that the MR can be located with as little as two independent TOA observations. 

104. The method of claim 103 wherein a method selected from the group 
consisting of the following methods is used to solve for the positional information of the 
MR: 

General Minimum Variance Unbiased Estimation, 
Best Linear Unbiased Estimation 
Maximum Likelihood Estimation, 
Least Squares Estimation, 
Method of Moments, 
General Bayesian Estimation, 
Linear Bayesian Estimation, and 
Kalman Filtering. 

105. The method of claim 103 where the TOA observations are based on 
estimating the RTD between the BST transmission and the reception of the MR's response 

30 to the BST transmission. 

106. The method of claim 89 in which Least Squares is used to solve for the 
positional information of the MR and where the number of BSTs used in locating the MR 
is optimized. 

35 



20 
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107. The method of claim 106 where ihe optimization of the number of BSs used 

in locating the MR is accomplished based on minimizing HDOP while maximizing the 
average received RSSI from all BSTs. 

5 108. The method of claim 89 in which Least Squares is used to solve for the 

positional information of the MR, the method comprising of a closed-form algorithm to 
provide an accurate initial position to start the Least Squares iterative process. 

109. The method of claim 108 where the closed-form algorithm is selected from 
10 the group consisting of LOCA, Plane Intersection, Bancrofi^s method, spherical 

interpolation, Schau and Robinson's method and Chan and He's method to provide an 
accurate initial position to start the Least Squares iterative process. 

1 10. The method of claim 88 in which the mean coordinates of the participating 
1 5 BSs are used as the initial position for Least Squares in the event that the closed-form 

causes Least Squares to diverge. 

in. The method of claim 88 in which the effect of geometry on AOA 

positioning is quantified using the design matrix in equation (45). 

20 

112. The method of claim 88 in which the effect of geometry on position 

estimation, with both AOA and TDOA observations, is quantified using the design matrix 
in equation (46). 

25 113. A method of estimating the location of a MR from a plurality of signals 

transmitted from a plurality of BSTs, which quantifies the effect of geometry on speed and 
direction of travel estimation using the design matrix in equation (51). 

114. The methods of estimating the location of a MR in claim 90 where the 

30 observations are selected from the group consisting of 
TOA, 
POA, and 
I- OA. 
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115. A system lo esiimaie ihc location of a MR from a plurality of signals 

transmitted from a plurality of BSTs, where misclosures and/or standardized residuals are 
used to flag observations that might contain a blunder. 

5 116. The system of claim 1 15 in which location information from each BST is 

weighted based on its corresponding RSSI and wherein a method selected from group 
consisting of the following methods is used to solve for the positional information of the 
MR: 

General Minimum Variance Unbiased Estimation, 
1 0 Best Linear Unbiased Estimation 

Maximum Likelihood Estimation, 

Least Squares Estimation, 

Method of Moments, 

General Bayesian Estimation, 
1 5 Linear Bayesian Estimation, and 

Kalman Filtering. 



117. The system of claim 115 which uses a method selected from the group 
20 consisting of Chaffee's method and LOCA to detect solution bifurcation. 

1 18. The system of claim 1 15 in which the hybrid TDOA positioning model in 
equation (42), which yields misclosures and residuals for the individual BSTs, is used. 

25 1 19. A system to estimate the location of a MR from a plurality of signals 

transmitted from a plurality of BSTs, where a combination of AOA and TDOA 
observations are used such that the MR can be located with as few as two AOA 
independent observations and one TDOA observation. 

30 120. A system to estimate the location of a MR from a plurality of signals 

transmitted from a plurality of BSTs, where a combination of AOA and TOA observations 
are used such that the MR can be located with as little as one AOA observation and one 
TOA observation. 
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121. A system lo estimate the location of a MR from a plurality of signals 

transmitted from a plurality of BSTs, where a combination of TOA observations are used 
such that the MR can be located with as little as two independent TOA observations. 

5 122. The system of claim 1 16 or 121 in which Least Squares is used lo solve for 

the positional information of the MR and where the number of BSTs used in locating the 
MR is optimized. 

123. The system of claim 1 16 or 121 in which Least Squares is used to solve for 
10 the positional information of the MR, the method comprising of a closed-form algorithm to 

provide an accurate initial position to start the Least Squares iterative process. 

124. The system of claim 115 or 116 in which the mean coordinates of the 
participating BSs are used as the initial position for Least Squares in the event that a 

15 closed-form algorithm selected from the group consisting of LOCA, Plane Intersection, 
Bancroft's, spherical interpolation, Schau and Robinson's method and Chan and Ho*s 
solution causes Least Squares to diverge. 

125. The system of claim 1 15 or 116 in which the effect of geometry on AO A 
20 positioning is quantified using the design matrix in equation (45). 

126. A system to estimate the location of a MR from a plurality of signals 
transmitted from a plurality of BSTs, which qtiantifies the effect of geometry on position 
estimation, with both AOA 2md TDOA observations, using the design matrix in equation 

25 (46). 

127. A system to estimate the location of a MR from a plurality of signals 
transmitted from a plurality of BSTs, which quantifies the effect of geometry on speed and 
direction of travel estimation using the design matrix in equation (5 1 ). 
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